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ABSTRACT 


This  analysis  is  concerned  with  the  use  of  radio- fre- 
quency (rf)  current  in  producing  lesions  in  the  brain.  The 
technique  consists  of  passing  the  current  between  an  active 
electrode  implanted  in  the  brain  and  an  indifferent  elec- 
trode located  on  the  skull.   The  resistivity  of  the  tissue 
causes  electrical  energy  to  be  converted  into  heat  thus 
raising  the  tissue  temperature.   A  problem  that  is  sometimes 
encountered  is  uncontrollable  boiling  of  the  tissue  near  the 
probe  tip.   A  theoretical  model  which  considers  the  combined 
effects  of  conduction,  blood  flow,  metabolism,  and  rf  heat- 
ing is  presented  to  show  the  relative  effects  of  each  of  the 
terms  on  the  temperature  field.   It  was  found  that  conduc- 
tive effects  and  rf  heating  effects  are  wiy    important  pa- 
rameters in  predicting  lesion  size,  while  the  blood  flow 
effect  is  only  marginally  important.   It  is  proposed  that  by 
maintaining  the  temperature  of  the  probe  at  a  low  constant 
value,  the  maximum  temperature  of  the  lesion  can  be  kept  well 
below  the  boiling  point,  yet  because  currents  and  time  can  be 
increased,  the  volume  of  tissue  destroyed  equals  or  even  sur- 
passes that  of  the  conventional  probe.   Finally,  the  radio- 
frequency  probes  are  compared  to  resistively  heated  probes. 
The  lesions  produced  by  these  probes  are  shown  to  be  signif- 
icantly smaller  than  the  lesions  produced  by  the  rf  probes. 
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I.   INTRODUCTION 

A.    STATEMENT  OF  PROBLEM 

In  order  to  treat  certain  disorders  occurring  in  the 
central  nervous  system  such  as  Parkinsonism,  unbearable 
pain,  and  malignant  tumors,  surgeons  at  times  find  it  neces- 
sary to  surgically  destroy  discrete  regions  of  nervous  tis- 
sue.  Numerous  methods  have  been  used  by  clinicians  to 
accomplish  tissue  destruction.   One  such  method,  which  has 
met  with  considerable  success,  employs  a  long,  thin,  needle- 
like probe  which  can  be  made  to  emit  a  low-power,  continuous 
wave  radio- frequency  current.   The  probe  itself  serves  as 
.the  emitting  electrode  and  a  second,  indifferent  electrode 
connected  to  some  remote,  external  point  on  the  body  acts 
as  a  ground.   As  current  is  passed  between  the  electrodes, 
the  resistivity  of  the  tissue  causes  electrical  energy  to 
be  converted  to  tissue  internal  energy  thus  giving  rise  to 
an  increase  in  local  tissue  temperature.   A  lesion,  or  re- 
gion of  destroyed  tissue,  is  formed  when  the  local  tissue 
temperature  is  raised  above  approximately  55°C  [3]. 

The  major  drawback  of  the  radio- frequency  probe  has 
been  excessive  heating  of  the  tissue  in  the  immediate  vicin- 
ity of  the  active  electrode  resulting  in  gas  formation,  car- 
bonization of  the  tissue  and  adherence  of  the  coagulum  to 
the  electrode.   Gas  formation  is  many  times  accompanied  by 
minor  "explosions"  so  that  the  method  often  reduces  to  a 
crude,  brutal  form  of  mechanical  tissue  destruction  [35]. 
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In  order  to  aid  surgeons  in  predicting  and  creating 
well  controlled  lesions,  it  is  necessary  to  model  the  heat 
transfer  process  occurring  in  the  region  surrounding  the 
active  electrode.   Several  attempts  have  previously  been 
made  at  formulating  such  a  model  [3,  35];  however,  solutions 
which  consider  the  combined  interaction  of  tissue  conductive 
effects,  blood  flow  effects,  metabolic  effects,  and  radio- 
frequency  heating  effects  on  the  change  of  tissue  tempera- 
ture with  time  have  not  yet  been  put  forth.   Such  a  model 
is  presented  in  this  thesis  for  spherical  and  cylindrical 
probes.   The  resulting  equations  are  solved  using  a  finite 
difference  technique  and  presented  in  a  form  that  can  be 
readily  used  by  surgeons. 

As  a  proposed  method  of  preventing  excessive  heating 
of  the  tissue,  a  constant  temperature  radio- frequency  probe 
and  a  high  resistance  heating  probe  are  also  analyzed  and 
discussed. 

An  additional  objective  is  to  present  a  summary  of  the 
various  techniques  that  are  used  to  create  lesions  in  tis- 
sue.  These  techniques  may  be  broken  into  several  broad 
categories:   mechanical,  chemical,  ultrasonic,  radio-active, 
x-ray,  freezing  and  electrical  methods  of  destruction.   A 
description  of  these  methods,  complete  with  the  advantages 
and  disadvantages  of  each,  follows. 
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B.    METHODS  OF  PRODUCING  LESIONS 

1 .  Mechanical  Destruction 

Perhaps  the  most  conventional  and  familiar  surgical 
technique  is  that  in  which  the  surgeon,  through  use  of  a 
scapel,  manually  makes  an  incision  and  removes  or  destroys 
unwanted  tissue.   This  technique  has  been  refined  for  use 
in  the  brain.   A  wire-like  knife  with  a  curved  or  angulated 
tip  is  concealed  in  a  suitable  length  of  hypodermic  needle 
tubing  of  approximately  lh   mm.  diameter.   The  needle  is  in- 
serted in  the  brain  to  the  required  depth,  and  the  knife  is 
pushed  out  and  rotated  360°,  thus  producing  a  lesion  in  the 
form  of  a  body  of  revolution.   Glees  and  others  [17]  have 
used  this  technique  and  note  that  the  lesions  produced  are 
well  localized  but  many  Limes  are  largo,  irregularly  shaped 
and  lack  uniformity.   In  addition,  Carpenter  and  Whittier, 
[6]  report  that  "The  lesions  appeared  to  be  dependent  upon 
interruption  of  local  vascular  patterns  and  to  a  variable 
extent  upon  the  blood  elements  liberated  into  the  tissue." 
As  a  result  of  this  dependence,  the  "rotating  knife"  method 
has  not  proven  successful  in  producing  consistent,  uniform 
lesions . 

More  obvious  disadvantages  of  this  mechanical  method 
include  profuse  hemorrhaging,  incidental  tissue  destruction, 
and  infection. 

2 .  Chemical  Destruction 

■Localized  necrosis  (tissue  destruction)  is  produced 
chemically  by  injecting  a  few  drops  of  zinc  chloride,  chromic 
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acid,  or  similar  corrosive  chemical  locally  into  the  brain. 
The  chemical  tends  to  diffuse  rapidly  and  large,  irregularly 
shaped  areas  of  tissue  destruction  are  produced.   Hyper- 
excitability,  convulsions,  and  unusual  pain  are  reported 
to  follow  lesions  produced  by  this  method.   In  addition, 
injected  solutions  tend  to  "leak"  out  around  the  needle  and 
cause  enlarged  needle  tracts  and  occasional  cortical  injury 
[6]. 

Because  of  the  unpredictability  of  the  size,  shape 
and  histological  characteristics  of  the  lesions  produced  by 
corrosive  chemicals,  this  method  is  of  little  value  for  use 
in  the  brain. 

3 .   Ultrasonic  Destruction 

Ultrasonic  waves  are  mechanical  vibr3tions  in 
solid,  liquid  or  gaseous  mediums  lying  above  the  range  of 
human  hearing.   The  piezo-electric  effect  is  adapted  to  the 
generation  of  ultrasonic  waves  at  frequencies  above  500 
kilocycles  (kc) ,  the  range  most  often  used  in  biological 
work.   Tissue  is  destroyed  primarily  by  microscopic  mechan- 
ical motion  although  heating  is  also  present.   Wood  and 
Loomis  [38]  give  an  account  of  the  effects  of  high  fre- 
quency sound  waves  upon  living  matter. 

The  ultrasonic  waves  are  produced  by  a  quartz 
plate  vibrating  in  resonance  with  a  high-frequency  current 
passing  through  the  plate.   It  is  possible  to  focus  these 
waves  at  a  central  focal  point  by  using  a  concave  quartz 
plate.   One  can  then  concentrate  as  much  as  150  times  as 
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much  ultrasonic  energy  at  the  focal  point  as  is  found  in  a 
region  near  the  vibrating  plate  [20].   As  a  result,  maximum 
tissue  destruction  occurs  at  the  focal  point  while  the  in- 
tervening tissue  remains  relatively  undamaged. 

Experiments  by  Lynn  and  others  [20]  show  that  the 
ultrasonic  technique  causes  an  area  of  destruction  on  the 
scalp,  followed  by  a  region  of  normal,  undamaged  tissue, 
and  finally  a  well  defined  focal  point  of  severe  destruction 
They  also  found  that  once  a  lesion  is  initiated,  it  appears 
to  extend  rapidly  and  progress  to  non- irradiated  parts  of 
the  brain. 

Presently,  the  successful  use  of  the  ultrasonic 
technique  is  hindered  by  the  scalp  injury.   Methods  are  now 
being  studied  to  eliminate  this  surface  destruction. 
4 .   Radioactive  Destruction 

Ey  implanting  a  small  glass  radon  seed  containing 
approximately  1  millicurie  (mc)  into  a  region  of  nervous 
tissue  in  the  brain,  one  can  create  a  sharp,  circumscribed 
area  of  destroyed  tissue  of  approximately  4  mm  in  diameter. 
Information  concerning  radon  dosage  values  for  the  produc- 
tion of  appropriate  neural  lesions  is  included  in  a  paper 
by  Borison  and  Wang  [4] . 

Stein  and  Peterson  [33]  report  that  the  lesions 
produced  by  radon  seeds  are  circular  and  discrete;  however, 
they  are  variable  in  size  and  volume  of  tissue  destroyed. 
They  postulated  that  "The  volumes  of  tissue  destroyed  were 
not  proportional  to  total  radiation  (millicuries-hours) , 
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probably  because  of  differences  in  radiosensitivity  of  dif- 
ferent parts  of  the  brain." 

Not  enough  data  were  taken  to  develop  conclusions, 
but  the  experimental  lesions  produced  were  not  lethal  and 
did  not  result  in  physiologic  disturbances.   The  lesions 
were  also  well  localized  and  consistent  in  configuration 
and  histologic  characteristics.   Perhaps  the  major  drawback 
of  this  technique  at  present  is  a  lack  of  radiosensitivity 
information  for  the  various  locations  in  the  brain.   A 
knowledge  of  this  property  is  needed  to  effectively  predict 
lesion  size. 

5 .   X-ray  Destruction 

By  stereotaxically  orienting  multiple  x-ray  beams 

tk^onrrK  +-  l-"  O  f  Vn  1  -»-f-  "!'"•  +Vir>^vi-if  i    r-   o  1    1    >r        ■»~\/'"\r-r*-iV\1J"v  "f"0  4-  r\  r*  1  1  O 

111      UUtijl  V      1    i    -  .-..  .t     J  \    ^  .<      1  .    _i         j  !        «.  (      ._  '  *_   J   »    v^       V..*  \.         *..      L.     <-     <^4 J-     _»-      ;  j   .-    -k.-'     w  ■     ..         J.     i    '     -       '■  L-   O  -i-    \-/     w    <~   ._> 

the  intersection  of  the  beams  at  a  predetermined  location 
in  the  brain.   This  will  result  in  tissue  destruction  at 
the  focal  point  in  much  the  same  manner  as  the  previously 
discussed  ultrasonic  technique. 

Carpenter  and  Whittier  [6]  point  out  "The  use  of 
x-ray  in  stereotaxic  technique  appears  promising,  but  such 
elaborate  equipment  in  neuroanatomic  laboratories  is  rare." 
In  order  to  use  this  method,  a  detailed  study  would  have  to 
be  under-taken  so  that  stereotaxic  orientation  by  mere  loca- 
tion on  the  skull  could  be  accomplished. 

The  advantage  of  this  technique  would  be  that  pre- 
dictable and  controlled  lesions  could  be  made  in  the  brain 
without  opening  the  skull.   Presently,  the  therapeutic  beams 
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available  for  this  method  produce  excessive  scatter,  and 
beams  small  enough  to  be  appropriate  in  an  experimental 
animal  are  not  of  sufficient  intensity.   There  are  x-ray- 
beams  that  would  produce  the  necessary  intensity  without  ex- 
cessive scatter,  but  their  use  for  stereotaxic  purposes  is 
unnecessarily  elegant. 
6 .   Cryosurgery 

Cryosurgery  has  proven  to  be  quite  useful  for  var- 
ious applications  in  the  brain,  especially  for  surgery  com- 
plicated by  bleeding  or  patient  discomfort  [15] . 

The  cryoprobe  is  a  hollow  tube  connected  to  a 
source  of  liquid  nitrogen.   The  probe  is  insulated  throughout 
.except  at  its  hemispherical  tip.   By  controlling  the  flow- 
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-196°C,  causing  the  tissue  in  the  immediate  vicinity  of  the 
tip  to  freeze  and  thus  be  destroyed.   The  cryosurgical  tech- 
nique is  bloodless  and  patient  discomfort  is  minimal.   In 
addition,  the  surgeon  can  produce  only  temporary  numbing  or 
can  destroy  a  region  of  tissue.   By  observing  a  patient's 
reactions  to  the  numbing  effect,  the  surgeon  can  ascertain 
whether  or  not  the  probe  is  in  the  proper  location  before 
lowering  the  temperature  to  the  lethal  level. 

A  description  of  the  lesion  formed  as  the  tissue 
freezes  is  given  by  Cooper  and  Trezck  [10,  12].   A  method  is 
also  presented  whereby  the  lesion  size  for  a  probe  of  arbi- 
trary dimensions  can  be  predicted.   The  rate  at  which  the 
temperature  field  develops  is  also  considered,  and  a  cryo- 
surgical aulas  is  proposed.   Both  steady  state  and  transient 
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situations  are  analyzed  in  their  work.   Another  heat  trans- 
fer model  of  the  cryoprobe  is  given  by  Rutkin  and  Barish 
[29]. 

Cooper,  Gionino,  and  Terry  [9]  describe  the  effects 
which  contribute  to  lesion  formation  by  freezing:   the  slow- 
ing down  of  the  blood  flow,  congested  blood  vessels,  dehydra- 
tion and  the  resulting  concentration  of  electrolytes, 
crystallization  and  the  consequent  rupture  of  cellular  mem- 
branes . 

Lack  of  hemorrhaging  during  the  operation  is  a 
marked  advantage  of  cryogenic  surgery.   The  small  amount  of 
hemorrhaging  that  does  occur  is  due  to  the  mechanical  pas- 
sage of  the  probe  through  the  brain  tissue.   However,  one 

large  probe  size  (2  mm.  or  greater). 
7 .   Direct  Current  (dc) 

It  is  possible  to  produce  generally  well  localized 
lesions  in  the  brain  by  applying  a  direct  current  between 
two  electrodes:   an  active  electrode  inserted  into  the  proper 
location  in  the  brain,  and  an  inactive  electrode,  or  ground, 
placed  on  some  other  position  in  the  body.   The  resistivity 
of  the  brain  tissue  causes  electrical  energy  to  be  converted 
to  heat,  thus  causing  the  local  tissue  temperature  to  be 
raised.   When  the  tissue  temperature  is  raised  above  55°C, 
the  tissue  is  destroyed. 

This  technique  can  be  used  with  the  active  electrode 
serving  as  either  a  cathode  or  anode.   When  the  electrode 
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acts  as  a  cathode,  the  lesions  formed  are  large  and  irregu- 
lar due  to  the  amount  of  gas  formed  at  the  tip  [6].   In  ad- 
dition, Mullan  and  others  [22]  found  the  cathodal  lesions 
to  be  hemorrhagic  and  hard  to  reproduce.   Similarly  the 
anodal  lesions  were  also  found  to  be  irregular  at  higher 
currents.   According  to  Mullan' s  work,  using  the  electrode 
as  an  anode,  surgeons  can  produce  small  and  discrete  lesions 
provided  they  slowly  increase  the  current  up  to  a  maximum 
of  2  milliamps  (ma)  or  prolong  the  duration  of  heating  up 
to  thirty  minutes . 

In  general,  however,  lesions  produced  by  the  direct 
current  method  vary  in  size,  shape  and  uniformity.   Watkins 
[37]  confirms  the  irregularity  of  the  d-c  lesions  and  adds 
that  the  lesions  sre  too  small  for  use  in  operations  or  man. 
He  also  states  that  "D-c  devices  suffer  from  an  electrode 
polarization  and  have  an  additional  disadvantage  that  the 
make- transients  and  break- transients  stimulate  muscle  and 
nervous  tissue."   The  polarization  of  the  electrode  is  an 
uncontrolled  variable  of  electrolysis  and,  along  with  the 
amount  of  gas  formed  at  the  electrode  tip,  is  a  cause  for 
the  lack  of  control  over  the  extent  of  tissue  destruction. 

The  fact  that  direct  current  stimulates  the  nerves 
is  an  advantage  in  one  respect:   the  surgeon  can  be  assured 
that  the  active  electrode  is  in  the  proper  location  in  the 
brain  by  observing  a  patient's  reactions  to  the  stimulating 
effect  of  the  direct  current. 

The  disadvantage  of  the  nerve  stimulation  is  more 
obvious.   Other  nerves,  as  well  as  those  directly  in  the 
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lesion,  could  be  stimulated  during  an  actual  operation. 
Furthermore,  should  the  current  find  its  way  to  the  heart, 
it  would  cause  the  heart  to  quiver  instead  of  pump.   In  ef- 
fect, the  surgeon  would  have  created  an  electric  chair  ef- 
fect . 

An  additional  disadvantage  of  the  dc  method  is  that 
direct  current  tends  to  have  preferred  paths  in  the  brain 
(paths  of  least  resistivity)  resulting  in  oval  shaped 
lesions  being  produced. 

One  final  disadvantage  of  the  direct  current  method 
is  noted  by  Herrero  [18]  :   "Animals  with  dc  lesions  usually 
had  extensive  scar  tissue  present  both  around  and  infiltrated 
.within  the  tissue."   This  scar  tissue  implies  that  the  blood 
vesseib  axOii^  wxtn  die  lisrvous  tissue  lias  uCCTi  \^.Z5 tr oyc»j. . 
Had  the  blood  vessels  remained,  the  body  would  more  easily 
accommodate  the  lesion. 

One  of  the  advantages  of  the  dc  method  over  most 
of  the  methods  discussed  previously  is  that  the  small  elec- 
trode can  be  painlessly  inserted  and  causes  minimal  damage 
to  the  surrounding  tissue  while  being  inserted  and  retracted. 
In  addition,  hemorrhaging  is  not  a  significant  problem  with 
the  dc  method. 

8 .   Radio-Frequency  (rf)  Current 

The  radio-frequency  technique  produces  lesions  in 
much  the  same  manner  as  the  direct  current  method.   An  ac- 
tive, needle-like  electrode  is  inserted  into  the  proper  lo- 
cation in  the  brain,  and  an  inactive  electrode  is  placed  on 
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the  skull  or  some  other  remote  portion  of  the  body  to  serve 
as  a  ground.   A  low-power,  continuous  wave,  radio- frequency 
current  is  passed  between  the  electrodes.   Tissue  destruc- 
tion is  accomplished  by  resistive  heating  of  the  tissue  it- 
self.  If  the  tissue  temperature  is  raised  above  55°C,  the 
tissue  is  destroyed.   The  factors  governing  the  amount  of 
heat  producsd  are  the  resistivity  of  the  tissue,  the  size 
of  the  electrode,  the  amount  of  current  used,  and  the  length 
of  time  the  current  is  applied.   An  inverse  relation  between 
current  and  time  gives  rise  to  the  amount  of  energy.   If  not 
enough  current  is  used,  no  tissue  destruction  occurs  for  an 
infinite  amount  of  time.   If  too  much  current  is  used,  vio- 
lent boiling  occurs  immediately  upon  switching  on  the  cur- 

Petty  and  others  [26]  conducted  a  study  correlating 
tissue  destruction  with  the  frequency  of  the  alternating 
current.   They  observed  that  low  frequency  current  resulted 
in  sporadic  or  no  tissue  destruction  at  all.   If  over  160 
ma  were  applied  at  four  (4)  kilohertz  (khz) ,  sparks  emanated 
from  the  electrode.   As  the  frequency  is  increased  above  45- 
khz  and  into  the  radio-frequency  range,  tissue  destruction 
proceeds  readily.   Rowland,  et.    at.     [28]  offer  an  explana- 
tion for  the  lack  of  coagulation  at  low  frequencies.   They 
suggest  that  coagulation  is  dependent  on  the  amount  of 
charge  per  pulse.   A  unit  of  charge  is  a  millicoulomb  (mcoul.) 
and  is  equivalent  to  a  milliamp  (ma)  times  a  second.   The 
total  energy  applied  by  an  alternating  current  is  determined 
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by  the  number  of  mcoul .  per  pulse  times  the  frequency  of 
pulses  per  second  times  the  length  of  time  current  is  ap- 
plied.  Recall  that  the  size  of  the  lesion  is  dependent  only 
on  the  current  (mcoul.  per  pulse  times  frequency)  and  the 
time.   Therefore,  a  low  frequency  current  applied  in  a  given 
amount  of  time  has  a  greater  charge  per  pulse  (mcoul.  per 
pulse)  than  does  an  equal  amount  of  high  frequency  current 
in  the  same  amount  of  time.   This  greater  charge  per  pulse 
could  possibly  be  what  causes  the  sparks  to  emanate  from 
the  electrode. 

In  reference  to  the  radio-frequency  technique, 
Aronow  [3]  explains  that  "For  surgical  purposes  it  is  nec- 
essary to  have:   (a)  precise  control  of  the  size  of  the 
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physical  position  of  the  lesion  in  a  relatively  inaccessible 
location;  and  (c)  precise  control  of  the  functional  position 
of  the  lesion." 

There  are  many  different  set-ups  of  the  radio- fre- 
quency apparatus  presented  in  the  literature  [1,  3,  13,  14, 
26,  31,  32,  34,  35,  36,  37].  Each  of  these  set-ups  contain 
certain  advantages  and  disadvantages. 

Excessive  heating  of  the  tissue  immediately  sur- 
rounding the  electrode  results  with  the  apparatus  proposed 
by  Aronow  [3] .   The  temperature  exceeds  100°C  and  the  tissue 
begins  to  boil  as  a  result.   In  addition,  there  is  gas  for- 
mation, carbonization  of  the  tissues  and  adherence  of  coagu- 
lum  to  the  electrode.   As  the  hard  coagulum  becomes  encrusted 
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on  the  electrode,  the  total  diameter  of  the  electrode  is  in- 
creased and  much  damage  to  the  surrounding  brain  tissue  re- 
sults as  the  probe  is  withdrawn  from  the  brain.   The  gas 
formations  are  actually  minor  "explosions"  so  that  the 
method  reduces  to  a  rather  brutal  mechanical  method.   Steam 
bubbles  are  formed  around  the  tip  of  the  electrode  and 
block  the  path  of  the  current.   An  audible  popping  sound  is 
an  indication  to  the  surgeon  that  boiling  is  taking  place. 
In  addition  to  the  obvious  disadvantages  that  occur  with 
boiling,  it  was  noted  that  patients  in  which  boiling  occur- 
red often  become  "somnolent  and  confused"  [35]. 

In  an  attempt  to  prevent  boiling  from  occurring, 
Van  den  Berg  and  Van  Manen  [35]  designed  their  "Coagrader" 
which  permits  controlled  and  predictable  lesions  to  be  formed 
"within  safe  limits  for  physiological  and  clinical  applica- 
tions."  In  addition  to  controlling  and  predicting  the  size 
of  the  lesion,  the  surgeon  must  be  certain  of  the  physical 
location  of  the  active  electrode  before  irreversible  damage 
is  initiated.   As  a  means  of  assuring  proper  probe  location, 
low- frequency  alternating  current  is  passed  between  the 
electrodes  and  the  patient's  response  is  analyzed.   Sias  and 
Goodman  [31]  point  out  that  low-frequency  alternating  current 
has  the  same  stimulating  effect  as  does  direct  current, 
whereas  radio- frequency  current  has  no  stimulatory  effect  at 
all. 

Another  means  in  which  a  surgeon  can  check  the  probe 
location  is  discussed  by  Brodkey,  et.    al .  [5].   They  discovered 
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that  a  temperature  range  of  40-49°C  can  reversibly  block 
nervous  function.   After  positioning  the  electrode  stereo- 
taxically  and  raising  the  temperature  of  the  tissue  sur- 
rounding the  electrode  to  the  reversible  range,  a  patient's 
response  will  indicate  whether  or  not  the  probe  is  directly 
on  target.   The  range  of  temperature  in  which  reversibility 
occurs  is  a  narrow  one,  however.   In  addition,  measurement 
of  the  temperature  of  the  tissue  surrounding  the  electrode 
is  not  a  simple  task.   It  is  easy  enough  to  find  the  tem- 
perature of  the  probe  tip  by  means  of  a  thermal  measuring 
device,  but  the  tip  temperature  is  not  necessarily  the  true 
tissue  temperature  surrounding  it.   It  was  shown  by  Dieckmann 
and  others  [13]  that  during  the  first  few  moments  of  heating, 
the  rate  of  increase  in  temperature  is  arpatpy  at  3  short 
distance  away  from  the  tip  than  at  the  tip  itself  due  to 
conductive  heat  losses  to  the  probe.   Even  externally  ap- 
plied thermocouples  have  been  reported  to  distort  the  tem- 
perature field  of  the  tissue  [5].   It  has  been  shown,  however, 
that  the  radio- frequency  technique  has  the  fine  control  nec- 
essary to  successfully  heat  reversibly  [5]  . 

One  problem  that  has-been  evident  with  the  r-f 
equipment  concerns  a  defect  in  probe  insulation.   Large,  ir- 
regularly shaped  lesions  have  been  known  to  form  at  points 
along  the  electrode  tract  due  to  current  "leaks."   Szekely 
and  others  [34]  developed  a  stylet  probe  to  minimize  this 
"leak." 

As  with  the  dc  probe,  the  very  small  size  of  the  rf 
probe  insures  that  the  overlying  brain  tissue  remains  undamaged. 
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If  the  probe  is  introduced  slowly,  vascular  structures  can 
be  pushed  aside  instead  of  punctured.   Closure  rapidly  takes 
place  as  the  probe  is  withdrawn,  and  discomfort  to  the 
patient  is  minimal  [39] . 

After  the  probe  is  inserted  in  the  brain,  and  the 
location  of  the  probe  has  been  checked,  appropriate  control 
of  the  physical  parameters  of  time,  tip  exposure,  and  tem- 
perature allow  well-circumscribed,  reproducible  lesions  to 
be  made  [36].   The  tip  diameter  also  affects  the  size  of 
the  lesion  produced. 

The  actual  stages  of  coagulation  are  discussed  by 
Petty  and  Edsall  [26].  There  are  three  stages  of  coagula- 
tion:  (1)  coalescense  around  the  electrode;  (2)  uniform 
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slowly  as  the  lesion  progresses);  and  (3)  boiling. 

Note  that,  unlike  the  dc  lesion,  the  rf  lesion  does 
not  follow  preferred  paths  in  the  brain  but  rather  unifor- 
mity in  all  directions  resulting  in  a  well -circumscribed 
spherical  lesion.   The  lesion  contains  three  zones:   a 
central  core  of  dense  dead  tissue  surrounded  by  a  zone  of 
damaged  tissue  and  finally,  a  narrow  zone  of  granulations 
which  clearly  mark  the  border  of  the  lesion  [2]. 

Mark,  et.    al .     [21]  point  out  that  the  blood  vessels 
running  through  a  lesion  are  often  times  spared.   "This  is 
partially  true  when  the  temperatures  used  at  the  tip  are  be- 
low 80°G."   Heating  above  100°C  results  in  rupturing  the 
blood  vessels.   As  was  noted  earlier,  this  corresponds  to  the 

boiling  temperature. 
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In  summary,  the  radio- frequency  method  allows  con- 
trolled and  predictable  spherical  lesions  to  be  produced 
within  safe  limits  for  physiological  and  clinical  applica- 
tions.  Hemorrhaging  is  not  a  problem  and  discomfort  to  the 
patient  is  minimal.   There  are  two  methods  available  with 
which  to  check  the  probe  location  before  irreversible  des- 
truction of  the  tissue  begins.   The  system  has  precise  con- 
trol over  the  volume  of  destruction  and  often  times  the 
blood  vessels  running  through  the  lesion  can  be  spared. 
Boiling  tends  to  be  a  problem  with  this  technique,  but 
proper  control  of  various  physical  parameters  can  avoid 
boiling . 
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II.   PHYSICAL  ANALYSIS 

A  region  of  nervous  tissue  in  the  brain  is  quite  com- 
plex.  It  contains  blood  vessels,  ventricles  and  is  nonuni- 
form in  composition.   The  grey  matter  in  the  brain  seems  to 
have  a  higher  resistance  to  temperature  than  does  the  white 
matter  [13].   Although  the  thermal  properties  (heat  capacity, 
thermal  conductivity,  density  and  thermal  diffusivity)  are 
known  to  within  an  accuracy  of  5%  [11] ,  there  are  still  many 
variables  in  the  brain  whose  exact  values  and  temperature 
coefficients  are  uncertain,  i.e.,  specific  resistance,  di- 
electric constant  and  circulation. 

Realizing  these  difficulties,  a  simplified,  one  dimen- 
sional heat  transfer  juudei  of  the  radio-frequency  probe  will 
now  be  developed.   Consider  the  brain  to  be  an  infinite, 
homogeneous,  imperfect  dielectric  of  resistivity  p'  and  di- 
electric constant  e.   The  active  electrode  is  considered  to 
be  a  conductive  sphere  of  radius  r0  embedded  in  the  brain, 
while  the  indifferent  electrode  is  a  spherical  conductor 
surrounding  the  brain  (see  Figure  1) . 

The  impedance  of  a  spherical  shell  of  material  of  width 
Ar  can  be  represented  by  an  incremental  capacitance  dC  in 
parallel  with  an  incremental  resistance  dR    [35]  (see  Figure 
2): 

dC   =    e    •    ^A      ;  dR   =   p.      Ar      • 


dZ    = 


Ar         '  K      4TTr2 

dR  dR 


1    +    jwdffdC        1    +    jcop'e 
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Figure  1.   Model  of  Spherical  Probe. 
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dC    = 

di?    = 
dZ    = 


£.  • 

4?rr2 

Ar 

P* 

Ar 
47rr2 

di? 

l  +  joodi?dC 


Figure  2.   Impedance  of  the  System 
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The  total  impedance  of  the  system  can  be  obtained  by  in- 
tegrating dZ  [35]  : 


oo  o 

r  dz  =    !_,_  r 

J  1  +  joop'e       J 


di?   = 


1  +  j  cop  '  e 


4iTr 


At  radio-frequencies,  however,  |jcop'ej<<l  for  brain  tissue 
[3,  35]  and  the  impedance  reduces  to  Z  =  p,/4iTr0.  For  ex- 
ample, for  brain  tissue: 


p'  -  150-200  ohm-cm   [35] 
e   *  70  e0   [35] 


sn  «  8.85X10"12   coulomb2   =  8#85xl()-14 


newt on -m2 


ohm- cm 


- 1 2    cpr 


Take 


p'  =  200  ohm -cm 


then 


p'e  =  12x10     sec 


Hence , 


joop'cj    =    12x10 


10 


CO 


where 


oo  is  in  Hz. 
Radio- frequencies  range  from  10  kHz  to  above  the  megahertz 
(MHz)  frequency.   Clearly,  even  in  the  megahertz  range, 
|jo)p'e|<<l,  i.e.  for  10MHz: 

oo  =  lOxlO6  Hz 
.  .  |joop'E;|  =  12x10"3<<1. 
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As  a  result,  in  the  radio-frequency  range  virtually  all  of 
the  current  passes  through  the  resistive  leg  of  the  circuit 
shown  in  Figure  2. 

To  solve  for  the  temperature  field  emanating  from  the 
radio-frequency  probe,  a  one-dimensional  energy  balance  is 
applied  to  a  spherical  element  of  tissue  (see  Figure  5) : 

Heat  into  element  by  conduction  -  Heat  out  of  element 
by  conduction  +  Heat  generation  due  to  metabolism  + 
Heat  input  due  to  blood  flow  +  Heat  generation  due   ^  ^ 
to  rf  heating  =  Increase  in  internal  energy  of  ele- 
ment . 

Equation  (1)  can  be  reduced  by  substituting  the  following 


L-V^  1  ilti> 


(i)  Cross-sectional    area   =   A   =    4Trr2 

(ii)         Elemental   volume    =    AV   =    47Tr2Ar 

(iii)      Heat    in   by    conduction    =    q.       =    -kA  -r—   =    -4irr2k   tt— 
v        J  }  lin  3r  dr 

da. 
in 
(iv)    Heat  out  by  conduction  =  q.   + 


dr 


2 


+  -I    d    ?in  Ar2 

2  !    3  r  2 

(Taylor  Series) 

(v)     Heat  generation  due  to  metabolism  = 

S  AV  =  S  4iTr2Ar 
°m      m 

(vi)         Heat    addition   due    to    blood    flow   = 
mbcb(Tb-T)AV    =   mbcb(Tb-T)4Trr2Ar 
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Volumetric  heat 
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Heating  due  to 

Ul-^J   -C  -i  -.  .  . 


m 


(by  conduction) 


Figure  3.   One-dimensional  Energy  Balance  on  Spherical 
Element . 
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(vii)    Heat  generation  due  to  rf  current  = 

47Tr2 
(viii)       Increase    in   internal    energy   = 

pc   |I  AV    =    pc   ~   47Tr2Ar. 

Note  that  in  the  model  the  effect  of  blood  flow  is  treated 
as  a  perfect  heat  exchanger;  that  is,  it  is  assumed  that 
the  blood  acquires  the  tissue  temperature  as  it  exits  a 
unit  volume  of  tissue.   This  is  consistent  with  a  model 
first  set  forth  by  Perl  [24,  25].   Rosomoff  and  Ho lad ay  [27] 
give  values  for  the  volumetric  blood  flowrate  m,  . 

If  the  metabolic  heat  generation  rate  is  assumed  to  be 

a        r*  r\T~\  c*  *f*  on  +•  t    +-        «~  o  t^        r\  /^        r*  om  r*  i   r\  r\  r\       t.t  i    +"  r\        *f-V»/~>        Kl/^/-\/^i         -T1    at.t       4-  o  "»~*^        *"*<? 

%cb(T0-T) 

where 

S 

rri       rp 


o     b   m,  c, 
b  b 

This  corresponds  to  redefining  the  reference  temperature. 
It  was  shown  by  Cooper  and  Trezek  []0]  that  in  the  case  of 
brain  tissue,  T0  is  on  the  order  of  0.5°C  higher  than  the 
systemic  arterial  blood  temperature  T,  .   For  the  case  of  a 
radio-frequency  probe,  the  blood  flow  convects  a  portion  of 
the  heat  generated  in  the  tissue  due  to  radio-frequency  ef- 
fects away  from  the  probe  tip. 

Upon  substituting  (i-viii)  into  (1) ,  the  following  equa 
tion  results: 
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3q.  T  2     i a 

-    -J^  Ar    +    \P   ;r   +   m,  c,  (Tfl    -    T)    4Trr2Ar 

(2) 

=    pc  |y  47Tr2Ar   +    0(Ar2)  . 

Dividing    (2)    by   Ar    and   taking   the    limit    as    Ar+O   yields 

-  ^f1  +  %&  +  mbVTo  -  T)  4"2  ■  p=  li  4"2-  ™ 

Note    that 

ST 

q.       =    -4-rrr2k   ^c—    • 


Hence , 


8r 


,  2«  i 

+ 


^-2^-  +  m,c,  (T0  -  T)47rr2  =  pc  ~   4ttt2  .  (4) 
4iTr2     b  bK    °     y  dt      v  J 


If  k  is  assumed  to  be  constant  and  (4)  is  divided  by  4-rrr2, 
the  equation  becomes 

Jr29l) 

\         \      9r/   +  J2^'-  +  m  c,  (T   -  T)  =  pc|£  .   (5) 

r2      3r       16TT2r'+     b  b  ^  o     '         H  9t     K    J 

It  is  now  convenient  to  introduce  the  following  non- 
dimensional  variables: 

m  c  r 
R  =  —  (Distance)   6  -   bt b     (Blood  Flow  Term) 
ro  K 


T     = 


-2-    (Time)  G    =    16    *rP^kT      ^rf  HeatinS  Term) 


ro  *  0  *"  0 


'  T  -  T 

o  B 

9  =  — fk (Temperature)  c}>  =  7=r 


Substituting  these  variables  into  (5)  yields 
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1 


R? 


3R 


dR 


66  + 


R" 


86 
8t 


(6) 


The  boundary  conditions  are 

(a)   at  the  probe  surface  (r=r0  or  R=l)  no  heat  is  con- 
ducted from  the  brain  tissue  back  into  the  probe,  or 


9_6 
9R 


R-l  =  °; 


(b)   at  large  distances  from  the  probe  the  temperature 
of  the  tissue  remains  constant;  that  is,  as  r->°° ,  T-*T0  ,  or 
in  non-dimensional  form,  as  R-*-»,  0->O. 

The  initial  condition  is  that  at  time  t=0,  T=T0  or  0=0 
for  all  locations  in  the  brain. 

Equation  (6)  along  with  the  boundary  conditions  and 
initial  condition  can  be  put  in  finite  difference  form  and 
solved  on  the  digital  computer.   For  the  details  of  this 
operation,  the  reader  is  referred  to  Appendix  A. 

A  similar  analysis  can  be  carried  out  for  a  cylindrical 

probe  resulting  in  the  following  equation: 

t 


R 


9_0_ 
8R/ 


R    dR 


36  + 


W 


22 

3x 


(7  J 


where 


IV 

47T2lzkT, 


0,  R,  and  3  are  the  same  as  above. 


The  boundary  conditions  and  the  initial  condition  are  the 
same  as  for  the  spherical  probe.   Equation  (7)  can  also  be 
solved  by  the  finite  difference  scheme.   The  reader  is  referred 
to  Appendix  B  for  the  details  of  the  cylindrical  solution. 
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Note  that  for  both  geometries  the  rf  heating  term  dies 
off  rapidly  as  the  radius  increases.   This  implies  that  the 
majority  of  the  rf  heating  is  generated  close  to  the  probe 
tip . 

As  was  mentioned  in  the  Introduction,  a  problem  that 
exists  with  conventional  rf  probes  is  tissue  boiling  near 
the  tip.   If  the  tissue  temperature  is  raised  above  approxi- 
mately 100°C,  the  rf  technique  reduces  to  no  more  than  a 
crude  mechanical  technique.   It  will  be  shown  that  even  for 
moderate  currents,  unwanted  boiling  may  occur.   In  order  to 
prevent  vaporization  at  the  tip  of  the  probe,  it  is  proposed 
that  a  constant  temperature  probe  be  developed  such  that  the 
maximum  temperatures  may  be  lowered  and  moved  away  from  the 

, I.    „        *-  -  „  TM, T   .  .  J-  1 J .        *-i.    •    ~        ~  1-   1 -     _         _1-j Z   .-    ~    1        _•--         u.1.   _ 

pXOUvJ         Ulp.  XllC        iUJ.ULi.UJl        L.U         tll-LO        JJiUL/XOIil       X  D        UUUalllCU        Xli        UilC 

same  manner  as  for  the  "insulated"  sphere  except  that  one 
boundary  condition  is  changed,  namely  the  temperature  at  the 
probe  tip  is  now  treated  as  a  constant.   The  results  of  this 
method  are  discussed  in  the  next  section. 

Another  technique  used  to  destroy  nervous  tissue  in  the 
brain  employs  a  "high-resistance"  probe.   In  this  method  the 
current  heats  the  probe  tip  rather  than  the  tissue,  resulting 
in  heat  diffusing  through  the  tissue  by  conduction.   This 
type  of  probe  is  discussed  by  Gildenberg  [16].   One  can 
solve  for  the  temperature  field  by  solving  equation  (6)  with 
the  rf  heating  term  equal  to  zero  and  using  the  constant  tem- 
perature boundary  condition.   Therefore,  the  equation  to  be 
solved  is 
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1   dj K  9R/    ftfi  _  89 

This  equation  can  also  be  easily  solved  using  the 
finite  difference  scheme.   However,  there  is  a  closed  form 
solution  to  this  type  of  problem  stated  in  Cars  law  and 
Jaeger  [7].   See  Appendix  C  for  details.   This  solution, 
when  compared  to  the  numerical  solution,  additionally  serves 
as  a  check  on  the  numerical  scheme  (Appendix  C) . 
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Ill .   RESULTS  AND  DISCUSSION 

In  a  previous  attempt  at  modeling  the  radio- frequency 
probe,  Aronow  [3]  assumed  that  "Distribution  of  heat  normally 
takes  place  primarily  by  convection  via  the  blood  flow, 
rather  than  simple  conduction."   He  then  postulated  a  model 
that  neglected  blood  flow  and  conduction  entirely.   It  was 
assumed  that  all  the  heat  generated  due  to  rf  current  went 
into  increasing  the  internal  energy  of  the  tissue.   The  sim- 
plified model  did  not  correspond  well  to  experimental  data. 
Aronow  indicated  that  neglecting  conduction  and  blood  flow 
are  the  most  probable  reasons  for  this  lack  of  agreement. 

Van  Den  Berg  and  Van  Manen  [3  5]  have  extended  Aronow' s 
model  somewhat.   They  include  the  conduction  term  along  with 
the  radio- frequency  term  but  neglect  the  blood  flow  term  in 
establishing  a  steady  state  solution.   They  conclude  that 
the  blood  flow  effect  is  negligible  and  that  the  model  is  a 
good  approximation  of  the  true  system.   Experiments  with  egg 
white  are  in  good  agreement  with  their  theoretical  analysis. 

As  was  seen  in  the  Physical  Analysis  Section,  the  pres- 
ent model  includes  the  combined  effects  of  rf  heating,  blood 
flow,  and  conduction  in  analyzing  the  transient  as  well  as 
the  steady  state  solutions.   Numerical  results  for  both  the 
spherical  and  cylindrical  probe  are  given  in  Tables  I-XVI. 
Values  for  6/G  are  given  at  various  non-dimensional  locations 
(R)  and. times  (x)  for  various  values  of  the  blood  flow  param- 
eter 3  (Tables  I -VIII).   Similarly,  values  for  8/0o  are  given 
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for  the  constant  temperature  cylinder  and  sphere  (Tables  IX- 
XVI)  . 

The  relative  effects  of  conduction,  blood  flow,  time, 
and  current  will  now  be  discussed.   In  addition,  the  constant 
temperature  probe  and  the  high-resistance  probe  will  be  com- 
pared to  the  insulated  probe  to  show  the  relative  advantages 
of  each. 

A.    THE  EFFECT  OF  CONDUCTION  OF  THE  TEMPERATURE  FIELD 

Figures  4  and  5  are  graphs  of  the  non-dimensional  tem- 
perature versus  the  radial   location  in  the  tissue  for  the 
blood  flow  parameter,  3,  equal  to  zero  and  non-dimensional 
time  t  equal  to  ten  (corresponds  to  a  real  time  of  approxi- 
mately 2.5  seconds  for  a  probe  with  a  radius  of  0.05  cm.). 
Two  curves  are  presented:   the  dashed  line  represents  Aronow ' s 
solution  which  neglects  conduction,  and  the  solid  line  repre- 
sents a  solution  which  includes  conductive  effects.   Note 
that  for  a  G  of  1.0  (corresponds  to  applying  a  current  of 
slightly  greater  than  2  ma  for  a  0.05  cm  probe),   0  is  10  at 
the  spherical  probe  tip  (R-1.0)  if  conduction  is  neglected 
(Aronow's  assumption).   This  value  is  well  in  excess  of  the 
boiling  temperature,  i.e.  boiling  occurs  at  T-100°C,  the 
initial  temperature  of  the  tissue  is  approximately  37°C, 
therefore 

a        =  Tboiling  "  T°    100-37  _  -,  7 
Uboiling       T0  37 

It  is  obvious  from  Figures  4  and  5  that  conduction  greatly 
reduces  the  temperature  at  the  probe  tip  compared  to  what  the 
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The  Effect  of  Conduction  on  a  Spherical  Probe, 
3  =  0.0,  t  =.10.0. 


41 


l&0fr 


6.*l 


0.2 


0.6 


—  Without  Conduction 

—  Conduction  Incluc/ed 


T    -    T 


2    r,    I 


IZP 


4'7T  2  1  2  kT 


T      = 


at 

V 


\ 


\ 


\ 


\ 


X 


\ 


J L 


|.o 


Figure    5 


_.    I         1 I L__1_J 

%Q      3.0      H.O       5.0      6-0  7-0       £.0      ?.£>      /0~O   M O 

R=  ^r0 
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3    =    0.0,    t    =    10.0. 
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tip  temperature  would  be  in  the  absence  of  conduction. 
Hence,  a  significant  portion  of  the  heat  generated  near  the 
probe  tip  is  rapidly  conducted  away  from  the  tip. 

B.    THE  EFFECT  OF  BLOOD  FLOW  ON  THE  TEMPERATURE  FIELD 

Figure  6  is  a  graph  of  the  temperature  response  in  the 
brain  due  to  a  radio- frequency  current  emitted  from  a  spher- 
ical probe  at  a  non-dimensional  time  of  ten.   Both  conductive 
effects  and  blood  flow  effects  are  included  in  this  graph. 
Recall  that  all  previous  models  assumed  the  blood  flow  ef- 
fect to  be  negligible.  'Note  from  the  graph  that  the  cooling 
effect  of  blood  flow  decreases  the  temperature  uniformly 
throughout  the  tissue.   As  the  blood  flow  parameter,  3,  is 
increased,  the  temperature  is  lowered  further.   Because  the 
blood  passing  through  the  lesion  has  a  lower  temperature  than 
the  heated  tissue,  heat  is  convected  out  of  this  region,  thus 
reducing  the  temperature.   To  gain  a  feeling  for  the  magni- 
tude of  this  effect  a  typical  value  of  3  will  be  calculated. 
A  typical  value  for  m,  c,  of  blood- flow  in  the  brain  is  0.01 
cal/cm3/sec/°C  [25]  and  a  typical  value  of  the  thermal  con- 
ductivity of  brain  tissue  is  0.0014  cal/cm/°C/sec  [11].   For 
brain  tissue,  therefore,  3  is  approximately  equal  to  7  r \ , 
where  r   is  the  probe  radius  expressed  in  centimeters.   A 
typical  radio- frequency  probe  might  have  a  radius  of  0.05  cm; 
therefore,  3  would  be  on  the  order  of  0.02.   As  can  be  seen 
from  Figure  6,  this  value  for  blood  flow  marginally  decreases 
the  temperature . 
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Figure  7  is  a  graph  of  the  blood  flow  effect  assuming  a 
cylindrical  probe.   Blood  flow  is  seen  to  play  a  slightly 
larger  role  in  Figure  7  than  Figure  6.   The  higher  values 
for  3  depicted  on  the  graphs  serve  to  illustrate  the  trend 
that  the  blood  flow  parameter  has  on  the  lesion  size.   Again, 
the  blood  flow  parameter  tends  to  decrease  the  temperature 
in  the  lesion  resulting  in  a  slightly  smaller  lesion. 

Figures  8  and  9  illustrate  another  effect  of  blood  flow. 
In  these  figures  the  temperature  at  the  probe  tip  (R=1.0) 
is  plotted  against  time  for  the  spherical  and  cylindrical 
probes.   At  low  times  the  blood  flow  effect  seems  to  be  neg- 
ligible, but  as  time  is  increased,  the  effect  of  blood  flow 
becomes  very  noticeable.   At  t  equal  to  10  it  can  be  seen 

tn  pt       *f-ho       effort       rv-P      l-\  1  o  /~\  ^       -Pi  m,r       i   c       *n  •"*  +       nDnl   i   rT-iKl   n  U  «~» -v*       o       v\-v**^l>^» 

of  radius  r0  equal  to  0.05  cm,  t  of  10  corresponds  to  approxi- 
mately 2.5  seconds.   As  the  time  is  prolonged,  the  blood  flow 
effect  is  increased  even  more.   The  larger  values  for  $  point 
out  how  increased  blood  flows  cause  steady  state  or  equilib- 
rium to  occur  at  an  earlier  time.   Again  blood  flow  is  shown 
to  decrease  the  size  of  the  lesion. 

Since  the  probe  radius  enters  the  blood  flow  term  as  a 
squared  term,  one  can  see  that  as  the  probe  radius  is  in- 
creased, the  blood  flow  effects  are  sharply  increased.   For 
a  larger  probe  the  blood  flow  effect  becomes  a  very  important 
parameter.   For  the  small  probes  and  small  times  used  in  the 
rf  technique,  however,  the  blood  flow  effects  are  only  mar- 
ginally important. 
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C.  EFFECT  OF  INCREASED  TIME  ON  TEMPERATURE  FIELD 

The  previous  section  pointed  out  the  fact  that  the  ef- 
fect of  blood  flow  was  increased  as  time  was  increased. 
Figure  10  indicates  the  trend  of  the  rf  temperature  field  as 
the  time  is  increased.   One  can  see  that  at  small  times,  the 
rate  of  temperature  increase  is  large,  but  as  times  increase, 
the  rate  of  increase  in  temperature  is  lessened  until  steady 
state  is  reached.   Steady  state  occurs  at  approximately  x 
equal  to  130  (corresponds  to  about  36  seconds  for  a  0.05  cm 
probe  radius) . 

D.  EFFECT  OF  rf  HEATING  TERM  ON  TEMPERATURE  FIELD  FOR 
INSULATED  PROBE 

Figure  11  illustrates  how  the  probe  tip  temperature  is 
greatly  increased  as  the  rf  current  term  is  increased.   Re- 
call that  a  value  of  0  of  1.7  corresponds  approximately  to 
the  boiling  point.   One  can  see  that  a  G  of  5  is  close  to 
this  critical  value  for  6  at  the  probe  tip  at  time  t=10.   If 
the  resistivity  of  brain  tissue  is  assumed  to  be  about  150ft 
cm  [35],  the  thermal  conductivity  cf  the  tissue  is  assumed 
to  be  0.0014  cal/cm/°C/sec  [11],  the  initial  temperature  of 
the  tissue  37°C  [3],  and  a  probe  of  radius  0.05  cm,  a  G  of 
5  corresponds  to  a  current  of  about  50  ma.   In  other  words, 
if  a  current  of  50  ma  is  applied  for  times  longer  than  about 
3  seconds,  boiling  will  occur.   In  order  to  prevent  boiling, 
a  much  lower  current  must  be  used,  (about  30  ma) . 

E.  CONSTANT  TEMPERATURE  PROBE 

It  is  proposed  that  boiling  can  be  avoided  alto. 
a  constant  temperature  probe  is  used  rather  than  an  insulated 
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probe.   Figure  12  is  a  graph  depicting  the  temperature  pro- 
file in  the  brain  for  various  values  of  the  radio-frequency 
heating  term  G  for  a  constant  temperature  spherical  probe 
(T-48°C,  3=0.1,  t=10.0).   As  G  is  increased,  the  maximum 
temperature  is  moved  away  from  the  probe  tip.   The  tempera- 
ture profile  is  more  flattened  than  it  was  with  the  insulated 
probe  (see  Figure  11).   Note  that  the  maximum  temperature, 
even  for  a  G  as  high  as  10,  is  well  below  the  boiling  point 
(8-1.7),  so  there  is  a  great  margin  of  safety  with  this 
probe.   Figure  13  gives  a  better  indication  of  this  safety 
factor.   In  this  graph  the  insulated  probe  is  compared  to 
the  constant  temperature  probe  for  a  given  value  3,  G,  and 
x.   The  insulated  probe  is  close  to  the  boiling  point  for  x 

the  boiling  point.   The  lesion  radius  depicted  by  the  in- 
sulated probe  for  t  of  10  is  approximately  0.16  cm,  i.e. 

TT   .    ^55°C,  T0  =  37°C 
Lesion       '   ° 

T        -  T 
Q        _   Lesion     °    55  -  37    n  An 

Lesion       T0  37 

From  Figure  13,  0-0 . 49n=>R=r/r0 -3 . 2 .   For  r0=0.05  cm,  r= 
0.16  cm  or  the  lesion  extends  (.16-. 05) cm  or  0.11  cm  from  the 
probe.   For  G=5.0,  $=0.1,  and  x=10,  the  size  of  the  lesion 
produced  by  the  constant  temperature  probe  is  considerably 
smaller.   However,  the  constant  temperature  probe  is  well 
below  the  boiling  point,  whereas  the  insulated  probe  is  dan- 
gerously close.   In  addition,  the  current  and  time  car  easily 
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be  increased  without  coming  close  to  the  boiling  tempera- 
ture.  For  example,  from  Figure  12,  it  is  evident  that  the 
constant  temperature  probe  (3=0.1,  t=10.0  and  G=10.0)  pro- 
duces a  lesion  of  approximately  the  same  size  as  one  de- 
picted in  Figure  11  for  the  insulated  probe  (3=0.1,  t=10.0 
and  G=5.0),  but  the  maximum  temperature  reached  by  the  con- 
stant temperature  probe  is  only  slightly  greater  than  70°C 
whereas  the  maximum  temperature  reached  by  the  insulated 
probe  is  over  92°C. 

Another  important  factor  can  be  observed  from  these 
figures.   Recall  that  if  the  tissue  temperature  is  kept 
below  80°C  (8-1.16),  only  the  nervous  tissue  in  the  lesion 
.is  destroyed  while  the  blood  vessels  running  through  the 
lesion  arc  spared  [21].   Note  that  for  all  values  cf  G  shewn 
for  the  constant  temperature  probe,  only  the  nervous  tissue 
is  destroyed,  whereas  for  G  greater  than  or  equal  5  for  the 
insulated  probe,  the  blood  vessels  as  well  as  the  nervous 
tissue  are  destroyed. 

One  final  point  to  note  in  regards  to  the  constant  tem- 
perature probe:   not  only  can  the  current  and  time  be  manip- 
ulated to  determine  the  size  and  maximum  temperature  of  the 
lesion,  but  the  temperature  at  which  the  probe  is  kept  can 
also  be  varied.   By  properly  combining  these  three  variables, 
optimum  results  can  be  accomplished. 

F.    HIGH  -  RESISTANCE  PROBE 

The  final  method  to  be  discussed  for  producing  lesions 
in  the  brain  is  the  high-resistance  or  pure  conduction  probe. 


5  5 


As  was  pointed  out  earlier,  this  technique  consists  of  pass- 
ing current  through  one  electrode  and  using  the  resistance 
of  the  electrode  to  cause  the  probe  to  increase  its  own  tem- 
perature, resulting  in  heat  being  conducted  away  from  the 
probe  into  the  surrounding  tissue.   This  heat  causes  a  tem- 
perature rise  in  the  tissue  and  eventually  tissue  destruc- 
tion results.   Figure  14  shows  the  resulting  temperature 
field  for  various  values  of  t  as  compared  to  the  insulated 
probe  temperature  field.   The  lesions  produced  by  the  high- 
resistance  probe  are  significantly  smaller  than  those  pro- 
duced by  the  rf  method,  but  it  should  be  pointed  out  that 
the  rf  lesions  have  an  added  amount  of  uncertainty  in  that 
the  electrical  properties  of  the  region  are  not  that  well 

VnnWll  l.'llprpac       •f-VtC»      Vi-inrVi--i~e»c-ic--f-o'nr-c»       lopinri       cinp       Jn«r-n   I  +       r!~- 

l\if  V    <      M.  X    y  .»,»»'_.    —     w     ■_*.  ^ .  Uliv  -  .     i      »_  -  j.  j  ^    -.-    u     L*  Ck  j  L  ■«--    ■ »  X   w  ■J    jl   \_.-  X  J.  .  ■    -i .    a_i   v„  vi   \*/  w   - .  >   li  L.  VX  "^ 

pend  at  all  on  these  electrical  properties. 

In  addition,  by  properly  controlling  the  temperature 
of  the  electrode  tip,  boiling  is  not  a  problem  with  the 
high-resistance  method.   In  summary,  this  method  produces 
safe,  predictable  and  very  small  lesions. 

G.    COMPARISON  OF  MODEL  WITH  EXPERIMENTAL  DATA 

An  attempt  was  made  to  check  the  proposed  model  with 
existing  experimental  data  [1,  3,  5,  6,  13,  28,  35,  36,  37]. 
However,  this  attempt  proved  unsuccessful  due  to  the  fact 
that  not  all  the  essential  parameters  were  presented  in  the 
reported  data.   In  order  to  correlate  the  analytical  results 
with  experimental  data,  it  is  necessary  to  know  the  size  of 
the  probe,  the  current  applied,  length  of  application, 
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temperatures  recorded  and  corresponding  positions  in  the 
brain.   At  least  one  of  these  parameters  was  either  missing 
or  stated  unclearly  in  the  reported  experimental  data. 

In  order  to  illustrate  this  fact,  a  sample  calculation 
will  now  be  given  from  the  data  collected  by  Aronow  [3]. 
The  data  consists  of  a  graph  of  temperature  rise  versus  time 
of  application  of  current,  and  a  table  which  gives  the  maxi- 
mum temperature  recorded,  the  distance  from  the  electrode 
at  which  this  maximum  temperature  was  recorded,  the  power 
used,  and  the  electrode  size.   However,  it  isn't  clearly 
stated  whether  the  electrode  size  refers  to  the  probe  radius 
or  diameter.   In  addition,  one  isn't  sure  of  the  current 
used  in  the  experiment.   Aronow  states  that  there  is  an  am- 

tabulate  the  values  of  current  or  voltage  used. 

It  will  be  assumed  that  the  electrode  size  refers  to 
the  probe  radius  therefore,  r0=0.05  cm.   A  temperature  rise 
of  45°C  is  reported  at  a  0.10  cm  distance  from  the  probe 
and  a  time  of  6  seconds.   Therefore,  the  following  data  is 
reported: 

r  0  =  0.05  cm 

r   =0.10  +0.05  =0.15  cm 

T   =  37  +  45  =  82°C 

t   =6  sec. 
It  will  be  assumed  that 

T0  =  37°C 

p'  =  200  ohm -cm 
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k  =  0.0014 


cm  C  s  e  c 


P  =  1  -g-3 

K      cm3 


c  =  °'9  g¥*C 

,  0.0014  — 6^al  „  2 

k  cm  Csec  n  nn^,    cm* 

a  =  —  = t-t:  -  0.0016  . 

PC  h  gm       _cal) 

I   cm3/     \gm  C/ 
Recall  the  following  non-dimensional  variables: 

r  -  r    b  -  mbCbr°2  x  -  at  e  -  T  "  T°  g  -  iy, 

An  attempt  will  now  be  made  to  substitute  the  previously 
mentioned  values  into  the  non-dimensional  variables  and  solve 
for  the  current. 

n      _      1  _      U  .1  J      CJU      _      7 

r0        0.05   cm 

T    -    T 

e  =  — _ — _  =  82"^7  a  i.22 


T0  37 


3    =   ^T''    .     I0'01    cT^sTc^05    ™,       ^    0>018 

0.0014 


cal 

Cs« 
cal 


2 


cmasec 


(0.OOI6   —  )  (6    sec 
at  \  sec/\  /  -z    oa 


r°  (.05cm) 


r  I2P'  SsT     •\/G    16ir2r    2kT0 

G   -    16Tr2r02kT0      :=5>I    =V I °_ 

pi 

Since  t  is  only  3.84  and  3  is  only  0.018,  the  effect  of 
blood  flow  will  be  neglected  in  this  example.   If  a  larger 
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probe  or  a  longer  time  were  used,  the  blood  flow  effect  would 
play  a  more  important  role. 

From  Table  I  for  an  insulated  sphere  at  R=3  and  t=3.84, 
it  can  be  found  that  6/G  0.078. 


G  =  __I_  =  L-|i  -  15.6 


Now, 


.078 


078 


i-V^ 


6ir  r   kT 


(15.6)  (16)  (tt2)  (.05cm)2  (.00 14 


cal 


cm  C  sec 


) (37°C) 


200  ftcm  • 


1  volt 

ampfi 


1  watt 


239  cal 


von  amp   watt  sec 


I  = V  67xl0"6 


amp ' 


I  = 


8.2  ma. 


This  value  for  the  current  appears  to  be  low  (Aronow 
states  that  a  typical  value  for  current  in  200  ma) .   The  ap- 
parent discrepancy  could  be  due  to  the  fact  that  the  probe 
size  is  in  fact  the  diameter,  or  the  distance  from  the  elec 
trode  is  considered  to  be  the  distance  from  the  center  of 
the  electrode  rather  than  from  the  outer  radius  of  the  elec 
trode.   If  the  current  had  been  reported,  and  the  vagueness 
in  the  above  terms  eliminated,  the  data  would  seem  to  be 
more  useful. 
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IV.   CONCLUSIONS  AND  RECOMMENDATIONS 

A.    CONCLUSIONS 

Surgeons  are  able  to  produce  controlled  and  predictable 
spherical  lesions  in  the  brain  with  little  or  no  hemorrhaging 
and  minimal  patient  discomfort  using  the  radio- frequency 
technique.   However,  care  must  be  taken  to  properly  control 
the  amount  of  current  and  length  of  time  the  current  is  ap- 
plied as  excessive  current  or  time  results  in  tissue  boiling 
and  uncontrolled  destruction. 

In  order  to  predict  lesion  size,  one  must  understand 
the  effects  of  conduction,  blood  flow,  metabolism,  and  radio- 
frequency  heating  on  the  tissue  temperature. 

Conduction  was  found  to  plav  a  verv  important  role  in 
determining  lesion  size.   The  effect  of  conduction  is  to 
distribute  the  heat  away  from  the  probe  tip  and  carry  the 
heat  into  the  brain.   The  conductive  effects  very  definitely 
cannot  be  neglected  in  attempting  to  predict  lesion  size. 

The  blood  flow  effect  was  found  to  be  marginally  im- 
portant in  determining  the  temperature  field  in  the  brain. 
Including  the  effect  of  blood  flow  tends  to  slightly  de- 
crease the  temperature  continuously  throughout  the  field. 
Thus,  the  blood  flow  convects  heat  away  from  the  region 
around  the  probe.   As  a  result,  steady  state  or  equilibrium 
is  reached  at  an  earlier  time  than  if  blood  flow  is  neglected 
The  effect  of  blood  flow  becomes  more  predominant  as  time  is 
increased.   Therefore,  if  heat  is  being  applied  for  any  sub- 
stantial amount  of  time,  the  effect  of  blood  flow  is  felt 
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re  strongly.   The  blood  flow  parameter,  3,  is  directly 
proportional  to  the  radius  squared.   As  a  result,  blood  flow 
plays  a  more  important  role  in  larger  probes  than  for  smaller 
probes.   For  the  probes,  currents  and  times  presently  used 
in  the  radio- frequency  technique,  the  blood  flow  effects  re- 
main marginally  important.   One  would  overestimate  the 
lesion  size  slightly  by  neglecting  these  effects. 

As  the  time  of  heating  is  increased,  the  lesion  size 
increases  at  a  decreasing  rate  until  steady  state  is  reached. 
For  a  spherical  probe  of  radius  0.05  cm  emitting  a  current 
of  about  20  ma,  steady  state  is  reached  in  about  40  seconds. 
If  a  larger  probe  or  a  greater  current  were  used,  steady 
state  times  would  be  larger. 

Increasing  the  current  also  results  in  higher  tempera- 
tures in  the  lesion,  especially  near  the  probe  tip  as  the 
majority  of  heating  occurs  in  this  region.   As  mentioned 
earlier,  excessive  current  levels  result  in  tissue  boiling, 
whereas  if  current  levels  are  too  low,  no  tissue  destruction 
occurs . 

If  the  temperature  at  the  probe  tip  could  be  maintained 
at  a  low  enough  constant  temperature  by  some  cooling  mecha- 
nism, both  the  current  and  the  time  can  be  significantly  in- 
creased, and  a  controlled  well- circumscribed ,  predictable 
lesion  would  result  without  the  ill  effects  of  boiling.  The 
maximum  temperature  may  be  moved  away  from  the  probe  tip  and 
kept  well  below  the  boiling  point  yet  high  enough  to  destroy 
nervous  tissue.   In  addition,  the  maximum  temperature  can  be 
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kept  low  enough  so  that  the  blood  vessels  running  through  a 
lesion  are  spared  while  the  diseased  nervous  tissue  is  still 
destroyed. 

The  final  probe  analyzed  in  the  thesis  is  the  high-re- 
sistance probe.   This  probe  is  also  kept  at  a  constant  tem- 
perature, but  the  current  is  used  to  heat  the  probe  rather 
than  the  tissue.   As  a  result,  the  conductive  effect  is  the 
only  mechanism  by  which  heat  is  introduced  into  the  surrounding 
tissue.   The  lesions  produced  by  this  method  are  signifi- 
cantly smaller  than  with  the  rf  method.   In  practice,  lesions 
produced  using  the  high  resistance  probe  should  be  easier  to 
predict  than  those  created  with  the  rf  probe  since  only  the 
thermal  properties  of  the  tissue  need  be  known. 

B .    RE  C OMM H N D AT 1 0 N  S 

It  is  recommended  that  experimental  work  be  carried  out 
with  both  radio- frequency  probes  and  high  resistance  probes 
to  check  the  validity  of  the  theoretical  models.   A  thorough 
literature  search  revealed  sparse  data  reported;  and  these 
data  proved  to  be  insufficient  to  substantiate  the  model. 
Not  enough  parameters  are  given  when  data  are  presented.   In 
order  to  effectively  compare  results,  it  is  necessary  to  in- 
clude the  following  parameters  when  recording  data:  size  of 
the  probe,  current  applied,  time  of  application,  temperatures 
recorded  and  corresponding  positions  in  the  brain. 

It  is  highly  recommended  that  the  constant  temperature 
radio- frequency  probe  be  studied  further.   In  particular, 
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simple  means  of  cooling  the  probe  tip  without  greatly  in- 
creasing the  size  of  the  tip  must  be  found.   This  probe 
would  appear  to  hold  great  promise  in  brain  surgery. 
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APPENDIX  A:   FINITE  DIFFERENCE  APPROXIMATION  FOR  A  SPHERICAL 

PROBE 

Statement  of  non-dimensional  equation  with  boundary- 
conditions  and  initial  condition: 


1_ 

R2 


3  R' 


8_0 

8R 


8R 


30  + 


R" 


s_e_ 


(A.l) 


(a) 


3J9 


R=l 


=  0. 


(b)  As  R+<»,  9+0. 

(c)  At  t=0,  6=0. 


where 


r         mbCbro 

R  =  —    3  =  -A-2 — 
r0  k 


x  = 


T  -  T 


G  = 


1-P' 


at 


* 


G  * 


There  seems  to  be  no  closed-form  solution  available  for 
this  equation  so  an  explicit,  central  difference,  finite  dif- 
ference technique  was  used  to  solve  for  the  temperature  field 
A  discussion  of  finite  difference  representation  of  various 
terms  is  given  by  Ozisik  [23] .   He  examines  the  finite  dif- 
ference approximation  of  the  terms  by  using  Taylor  Series 
expansions . 

Consider  a  spherical  coordinate  system  with  a  region 
interlaced  with  a  R,t  network  such  that  the  coordinates  of 
a  nodal  point  (i,n)  is  given  as  R=iAR  (non-dimensional  radial 
location)  and  x=nAi  (non-dimensional  time) .   The  non-dimen- 
sional temperature  at  time  t=uAt  and  radial  location  R 


defined  as  G 


n 
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Referencing  Ozisik  [23]  the  following  transformations 

can  be  made : 

(A. 2 

V20  =  J- 


R 


9R 


,n 


36  =  36.  . 
l 


AR2 


-ikUY (^U)\- ^l 


(A.  3) 


G 

ii 

9t 


(iAR)  k' 

en+1  -  en 
i i 

At 


(A. 4) 


(A.  5) 


Substituting  these  finite  difference  approximations  into 
(A. 1)  yields  • 


fc  (Bk-!  *  C1+r>?+i} 


AR 


en 


e:       -   ev 
i_ i 

Ax 


Solving    for    6 


n+1 


6n+1 

l 


At 


AR' 


6-0°ti +  Hk+i 


B6^   + 


_G 

(iAR)1* 


+  u 


2At 
AR2 


(A. 6) 


(A. 7) 


^ei  +  T^ 


Let   M   = 


At 
AR2 


(A.  8) 


^=M(H>i-l+(1+l)il)    +     {l-2M-3AR2M 


>Qi    +    TiARj^ 


The  stability  criteria  serve  as  a  constraint  on  the  prob- 
lem.  A  stability  condition  must  be  introduced  to  limit  the 
maximum  size  of  time  steps  that  can  be  chosen  for  a  fixed  AR. 
The  stability  criteria  for  this  problem  will  be  stated  without 
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proof:   1  -  2M  -  3AR2M  must  be  >  0  (the  interested  reader 
is  referred  to  Ozisik  [23]). 


then 


If  1  -  2M  -  3AR2M  =  0 


M  = 


2  +  3AR2  ' 

Recall  that  M  =  At/AR2 .   Choosing  AR  automatically  specifies 
At,  i.e.  Ai  =  MAR2.   Hence, 


GMAR2 


Recall  from  (A.l)  that  the  boundary  conditions  are 


(A.  9) 


O)  30 
3R 


R  =  l 


=  0 


(b)  As  R+00,  0+0. 

and  the  initial  condition  is 

(c)  At  t=0,  6=0. 

Using  the  finite  difference  scheme  the  boundary  conditions 
become 


fin     -fln 
6i*  +  l   6i»-l 

2AR 


0   or   9.*  i  =  8.*.! 

i;,-l     i  *  + 1 


where  i*  corresponds  to  R=l . 0=i*AR , 


(A. 10) 


)n 
JN+1 


6. 


n 


N -1    n      nn 
2AR =  °   or   GN-1 


where  i  =  i*,  i*+l,  ...  ,N-1,  N. 
The  initial  condition  becomes 


6.  =  0.0. 

i 


=  6 


n 

N+l 


(A. 11) 


(A. 12) 
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Summary : 
Differential  equation 


^+1=M{(l-i)e?.1  +  (l4>?+1)^^-         CA-13J 

Boundary  conditions 

(a)  en    =  Qn 
ui*-l    °i*+l 

(b)  Qn    =   n 

N-l     DN+1  * 

Initial  condition 

(c)  ec.  -  o.o. 
i 

Comput:er  Program  1  was  used  to  solve  this  problem. 
Richardson's  Technique  [8]  was  used  to  test  the  accuracy  of 
the  solution.   Three  numerical  solutions  usina  different 
sized  increments  of  AR  were  obtained.   Because  a  central  dif- 
ference approximation  has  a  discretization  error  of  the  order 
of  AR2  [23],  one  can  plot  the  approximation  to  the  tempera- 
ture versus  the  increment  squared  (AR2)  for  each  of  the. throe 
increments  used  and  extrapolate  the  plot  to  give  the  exact 
solution.   Figure  A.l  shows  Richardson's  technique  as  applied 
to  the  spherical  program.   An  increment  AR  =  0.1  proved  to  be 
a  convenient  interval  for  the  program  to  work  with  and  was 
well  within  the  limit  of  accuracy  desired  by  this  approximate 
solution . 

An  additional  problem  that  was  analyzed  for  this  thesis 
was  the  case  of  a  constant  temperature  probe.   Recall  that 
this  probe  serves  to  lower  the  maximum  temperature  incurred 
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at  the  probe  surface  and  thereby  serves  to  prevent  the  tis- 
sue from  boiling.   The  differential  equation  and  initial 
condition  remain  the  same  as  stated  in  (A. 13).   Boundary 
condition  (b)  also  remains  unchanged,  but  boundary  condition 
(a)  is  no  longer  valid.   The  probe  is  now  considered  to  have 
a  constant  temperature  rather  than  being  thermally  insulated, 
or 

B1}*   =  Constant.  (A.  14) 

Computer  Program  2  was  used  to  solve  this  problem. 
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APPENDIX  B:   COMPLETE  SOLUTION  FOR  THE  CYLINDRICAL  PROBE 

Statement  of  Problem: 

A  cylindrical  probe  of  radius  "r  "  is  embedded  in  tissue 
which  is  initially  at  a  uniform  temperature  of  T  .   Within 
the  tissue  there  is  volumetric  heat  generation  due  to  metab- 
olism of  strength  S  .   In  addition,  there  is  blood  flowing 

°  m  b 

through  the  medium  with  a  flowrate  nw  ,  specific  heat  c,  ,  and 
temperature  T,  .   Assume  that  the  length-to-diameter  ratio  of 
the  probe  is  large  enough  so  that  edge  effects  may  be  ne- 
glected and  assume  the  tissue  is  homogeneous  and  infinite  in 
extent.   A  radio- frequency  current  I  is  emitted  from  the 
.probe,  and  the  tissue  is  heated  due  to  the  resistivity  p'  of 
the  tissue.   Assume  that  the  heat  is  entirely  dissipated  in 
the  tissue  and  that  none  of  the  heat  is  conducted  back 
through  the  probe,  i.e.  the  probe  is  thermally  insulated. 
Also,  assume  that  at  an  infinite  distance  from  the  probe,  no 
heat  is  conducted.   Solve  for  the  probe  temperature  assuming 
constant  properties  and  one-dimensional  radial  heat  flow. 

Notation: 

Symbol  Description  Typical  Units 

a  Tissue  thermal  diffusivity  cm2/sec 

c  Tissue  specific  heat  cal/gm/°C 

c.  Blood  specific  heat  cal/gm/°C 

dR  Incremental  tissue  resistance  ohm 

I  Current  amps 

k  Thermal  conductivity  of  tissue  cal/cm/sec/c 
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Symbol  Description  Typical  Units 

m.  Volumetric  blood  flowrate  gm/cm3/sec 

p1  Resistivity  of  tissue  ohm-cm 

p  Density  of  tissue  gm/cm3 

r  Radial  location  in  brain  cm 

r  Probe  radius  cm 

S  Volumetric  metabolism  rate  cal/cm3/sec 

111 

t  Time  sec 

T  Tissue  temperature  °C 

T,  Temperature  of  blood  °C 

T  Initial  tissue  temperature  °C 

Non-dimensional  groups: 

Temperature         Distance  Time 

rp        rp 

1      x  0  -r  nt 

6 ^ R  -  ~  T  -  ~ 

1  0  0  o 

Blood  flowrate  rf  heating  term 

ft    -   mbCbr°2  •  f,    _         IV 

p  k  b  4TT2l2kT0 

Make  a  one-dimensional  energy  balance  on  an  element  of  tis- 
sue (see  Figure  B.l):   He'at  into  element  by  conduction  - 
Heat  out  of  element  by  conduction  +  Volumetric  heat  generation 
due  to  metabolism  +  Heat  input  due  to  blood  flow  +  heat  gen- 
eration due  to  rf  current  =  Increase  in  internal  energy  of 
element . 

ST"  ?)  T 

Heat  in  by  conduction  =  q.   =  -kA^r-  =  -k(rA6AZ)  tt— 
7  lm      3r      v     J  3r 

Heat  out  by  conduction  =  q.   +  —  ~-    Ar  +  -~-r      „■  2 —  Ar2+  .... 
}  nin     9r        2 !   dr^ 


(Taylor  Series) 
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Volumetric  heat  generation  due  to  metabolism  = 

S  AV  =  S  (rA6ArAZ) 
m      m^       J 

Heat  input  due  to  blood  flow  = 

mbcb(Tb-T)AV  =  mbcb(rA6ArAZ)(Tb-T) 
Heat  generation  due  to  rf  heating  = 


I2d/?  =  I 


p  '  Ar 


rAGAZ 
Increase  in  internal  energy  of  element  = 

pc  ^  (rAGArAZ) 

O  L 

Combining  the  above 

9°i  •  T  2   1  A 

-  -^  A^  +   AQA7-  +  m,c,  (T,  -T)rA0ArAZ  +  S  (rAGArAZ) 
8r        rAGAZ      b  bv  b  J  mv        J 

_„~„aOa„a7.A('a^2'\  f"R   "I  "1 

or  by  dividing  by  Ar  and  taking  the  limit  as  Ar->0 


(B.2) 

■^   +  -t^tt  +  m,c,  (T,  -T)rA6AZ  +  S  rAGAZ  =  pc  -^  rAGAZ. 
dr  rAOAZ     b  bK    b  J  m        ^   3t 


+  — l-J^L-   +    m.  c.  fT.  -TlrABAZ  +  S  rAOAZ  =  nr    £1 


Recall  q.   =  -krAGAZ  j^- 

lin  3r 

3(krAGAZ  —)  2  , 

5 —  +  aaav  +  m,  cK(T,-T)  rAGAZ  +  S  rAGAZ 

8r  rAGAZ     b  b v  b  *                       m 

=  pc  ~  rAGAZ.  (B.3) 

a  X. 


Dividing  through  by  rAGAZ 


S   +   ?  aLa72  +  mi  ct  (TT  -T)  =  pc  |X.   (B.4) 
r  •   3r       m    r2AG2AZ2     b  b v  b   ^       9t 
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Heat  generation 
due  to  rf 
current 


q   . (by  conduction) 
nout^  '  J 


Volumetric  heat 
A2  /"generation  due 

sm 


'  /^generation  i 
J     to  metaboli 


Heating  due  to 
blood  flow 


q.  (by  conduction) 


Figure  B.l.   One-dimensional  Energy  Balance  on  a  Cylindrical 

Element . 
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For  a  constant  k  (B.4)  reduces  to 


k  Hr  dr  )         c       I2p'     m    „  fr      r,         rt„  8T 

—  * +  S    +    a  .  Q?  ?„  9   +  ITl,  C,  (Tv-T)   =  PC  rr  . 

r    8r        m   r2A62AZ2     b  b  *■  b   '       8t 


At  this  point  a  further  constraint  will  be  placed  on  the  sys- 
tem.  Prior  to  the  insertion  of  the  probe  it  will  be  assumed 
that  the  temperature  of  the  medium  is  constant.   In  addition, 
it  will  be  assumed  that  metabolism  and  blood  flow  have  been 
active  long  before  the  insertion  of  the  probe. 
If  this  be  the  case,  then 

S   =  m,c,  (T.-T) 
m    b  b v  °  } 

or 

o     b    m'  b  b 

This  implies  that  the  medium  temperature  is  warmer  than  the 

fluid  by  an  amount  S  /m,  c,  prior  to  the  insertion  of  the 
;  m   b  b  r 

probe.   At  steady  state  the  fluid  thus  removes  the  heat  added 
by  the  volumetric  heat  source. 

Putting  Equation  (B.5)  in  the  form 

k  Hr  3r/t      (L  Sm  ^   T\  x    I2p'       „  8T    ,R  ,, 

^ +  mi  ci  ?  \Ti  -  -  T  (  +   ,.n,A7-7  =  pc  ttt-  .  (B .  6) 

r    3r      b  b/\b    m,  c,  /.    ;    r2AG2AZ2       8t       J 

suggests  substituting  T\  for  (T,  +  S  /m,  cu) .   When  this  is 
&&  °   o      *•  b    m  b  b' 

done,  Equation  (B.6)  reads 

dfv   —) 

IS.        ^      8r  '      +   m    r    TT    -Tl    +  l2P*       =    oc   —  (B    71 

r  9r  mbCbUo    U         r2A92AZ2         pc    dt    '         l        } 

Let   A6    =    2tt    and   AZ    =    1    to   get 
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k  3(r  It) 


dr 


*   mbcb(T0-T)    +    .42irgI2      =PCH.  (B.8) 


Now,    to    cast    (B.8)    into   non-dimensional    form   let 

R  -  I_.      b  =  mbCb"f°     x  =  SLL.  e  =  T  I  T°     g'   -       i2p' 

R        r0         3  k  T         r/  To  4u2l^kT0 

Substituting  the  above  non-dimensional  quantities  into  Equa- 
tion (B . 8)  gives 

1  HK_3_RJ/_        G«  _  8G  f    . 

Equation  (B.9)  is  the  equation  that  must  be  solved  in 
order  to  describe  the  temperature  of  the  tissue  as  a  func- 
tion of  time  when  a  radio-frequency  current  is  emitted  from 
a  probe  embedded  in  a  medium  with  a  volumetric  heat  source 
and  a  flowing   fluid.   In  order  to  solve  the  problem,  bound- 
ary conditions  and  an  initial  condition  must  be  stated.   The 
boundary  and  initial  conditions  are  as  follows: 

(a)  None  of  the  heat  generated  due  to  rf  current  in 
the  tissue  is  conducted  into  the  probe  (probe  is 
thermally  insulated) ; 

AT  r  =  r0 ,  f£  =  0 

or,  in  non-dimensional  form 

at  R  =  1,  99/8R  =  0; 

(b)  as  r-*00 ,  the  medium  remains  unaffected 

as  r-*»   T->T 
o 

or,    in   non-dimensional    form 
as    R->°°  ,    9-+0  ; 
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(c)   at  time  t  =  0,  T  =  TQ 

or,  in  non-dimensional  form 
at  t  =  0,  9  =  0 . 
Summary 

Differential  equation 


i  »(»  ft) 


G'    36 


R    3R      "  B9  +  RT=  3T  •  (B-10> 
Boundary  conditions 

'  9rIr=1  =  °  (B.ll) 

as  R-*»,  e+o.  (B.12) 


Initial  condition 


at  t=  0,  e  =  0. 


Solution 


An  explicit,  central  difference,  finite  difference  ap- 
proximation was  used  to  solve  for  the  temperature  field. 
Consider  a  cylindrical  coordinate  system  interlaced  with  a 
R,t  network  such  that  the  non-dimensional  temperature  at  a 

time  T=nx  and  radial  location  R=iAR  is  defined  as  6 . ,  where 

l ' 

n=l,2,...,  and  i=i*,  i*+l , . . . ,N-1 ,  N,  and  i*  is  the  location 
where  R=1.0. 

Referring  to  Ozisik  [23]  the  following  transformations 
can  be  made : 

^  -   J  4^-  ■  zir  {M9?-!  *KJe?*i  ■»;}»■") 
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36  =  36 


n 


(B.15) 


91 
R2 

M 

8t 


G' 
(iAR)2 

0n+1  - 

l 


0 


n 


At 


(B.16) 


(B.17) 


Substituting  these  finite  difference  approximations  into 
(B.10)  gives 


AR: 


{(i-zrki  *  fi^KiJ 


0n   + 


AR2      i 


G' 

(iAR)"2 


0n+1      -    en 
l i_ 

At 


(B.18) 


c  i  •    -r   An+1 
Solving  for  6-    gives 


n+1  At 


AR2 


{(.      lln 


/,  .    lLn       ( 


(B.19) 


f !       2/AT 

^1+2TJ°i+lJ  *  \  1_ART 


I  .n 


V 


G  :  At 

TTaX>2 


Let  M  =  At/AR2 


(B.20) 


'i+1-M{Me;-i  +  (1+^ei+i 


+  /  i-2m-3ar2m  re 


n  +  G'At 


i    (iAR)2 


To  maintain  the  stability  criteria  [23], 


then 


1  -  2M  -  3AR2M  =  0 


M  = 


2  +  3AR2 

Recall  M  =  At/AR.   Choosing  AR  automatically  specifies  At 
i.e.  At  =  MAR2.   Hence, 


(B.21) 


(B.22) 


(B.23) 


^■"fM'S-i^HHi^ 
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Recall  the  boundary  conditions  from  (B.ll)  and  (B.12) 
(a)  36 


3RlR=l  ° 

-  (b)  as  R-*»,  9+0. 

Using  the  finite  difference  scheme,  the  boundary  conditions 
become 


An  A11 

i*+l  i*+l  n  Mn  ~n 


2AR 


0  or      eV*^    =    9iA  +  1  (B.24) 


en      -  en 

N+l          N-l  n  nn  nn  rT>    or. 
2AR =    °           or      °N-1      =    Vl  ^B-25^ 

where  i  =  i* ,  i*+l,...,  N  -  1,  N,  and  i*  corresponds  to 
R  =  1.0.   The  initial  condition  becomes 


l 
Summary 
Differential  equation 


6.  =  0.0.  [B.26) 


•T1  -  «  {6-itRi  *  (-ir>i4  +  G-i 


(B.27) 


Boundary  conditions 

(a)  on,  i  -  en,+1 

(b)   n    _   n 

°N-1     dN+l  * 

Initial  condition 

(C).   6°  =  0.0. 

1 

By ■ specifying  $,  the  blood  flow  factor,  and  AR,  the  in- 
crement, M  can  be  found  and  Ax  can  be  found.   Computer  Progi 
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3  does  this  and  also  solves  for  the  temperature  at  various 
times  and  locations.   Richardson's  technique  [8]  was  used  to 
check  the  accuracy  of  the  program.   Figure  B.2  shows  Rich- 
ardson's technique  as  it  applied  to  the  cylindrical  program. 

As  with  the  spherical  probe,  the  constant  temperature 
boundary  condition  was  analyzed  with  the  cylindrical  probe. 
The  only  difference  to  the  problem  as  stated  in  the  previous 
summary  is  boundary  condition  (a) .   This  boundary  condition 
becomes 

6?*  =  Constant.  (B.28) 

Computer  Program  4  was  used  to  solve  this  problem. 
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APPENDIX  C:   SOLUTION  FOR  THE  HIGH- RESISTANCE  SPHERICAL 

HEATING  PROBE 

Differential  equation 

i  d(RZ   FrJ    Rfi  _  dB 

Boundary  conditions 

At  R  =  1 ,  0  =  90  (a  constant)  (C.2) 

As  R-k»,  0*0.  (C3) 

Initial  condition 

4 

At  t  =  0,  6  =  0.  (C.4) 

Expanding  (C.l)  gives 

2   91  .  l!l  .  Rfi  =  3_i  rr  «-, 

XX     V  1\  U  i\-  0  t 

Using  a  change  of  variables  technique  as  described  in  Carslaw 
and  Jaeger  [7],  let  6  =  v/R. 
Then, 


(I) 


dQ_        3  \RJ  v_       Sv 

3R     8R        R2    R   3R*  l    J 


ili  =2v        _J_  9v        _1    3v    +    I  llv  rr    7^ 

9R2         R3         R2    8R        R2    8R        R    3R2     '  L    * /J 


or 


826    2v    2_  9v    1  82v 
3R2    IT3-   R2  8R    R  8R2  * 

Substituting  these  new  variables  into  (C.5)  yields 

2v    _2   9v    2v     2  9v  ,  1  32v    qX  =  1_  9v 
R1"    R2   SR    R3  '  R2  3R    R  3R2     R 

132vv_13v  fr    81 

R3RT"^    R    "    R    3?    '  l         J 
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Multiplying  through  by  R 

From  Car  si  aw  and  Jaeger  [7],  p.  134., 
equation  is 


(C9) 


the  solution  to  this 


v  = 


1  Voe-/B   (R-D  Lfc  IRiii  .  vpA  +  e/B"  (R-D 

'         2  yr  ) 


erfc 


(R-l) 
2-/t 


+  V&T 


(CIO) 


or 


6 


_1 
2R 


-  "VE"  (R-l)j 
e  v        Merfc 


(R- 

2 


^T 


+   eVTCR-lD    erfcfi^l     +  VBT  j 


Ceil) 


This  equation  was  solved  using  Computer  Program  5.   The  re 
suits  are  shown  on  Figure  14. 

The  steady  state  solution  can  be  obtained  directly: 


d2v 


wr  -    3v  =  0  ^t>  v  =  C,  e 


■^e  CR"1)  +  c  e~^  (R_1) 

L-i 


Applying   boundary    condition    (C.3) 
As    R+oo,    e->0  =J>C7    =    0 


•         • 


fl    =   _I  e'^   CR-1) 

6  R   e 


At  R  =  l,       e  =  e0  =^>c1  = 
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•    • 


_e  £^   (R"13 

e0  "  r 

A  comparison  of  the  high- resistance  program  and  the  finite 
difference  program  for  the  constant  temperature  probe  (tak- 
ing G  =  0)  serves  as  an  additional  check  on  the  finite  dif- 
ference program.   Figure  C.l  shows  this  comparison.   As  the 
graph  depicts,  the  two  solutions  fall  exactly  on  one  another, 
thus  giving  additional  assurance  to  the  previous  programs. 
When  comparing  the  lesion  size  of  the  high- resistance 
technique  with  the  lesion  size  of  the  insulated  rf  technique, 
it  is  desirable  to  compare  the  size  at  steady  state  values. 
In  order  to  do  this,  the  rf  probe  must  be  fixed  so  that  it 
has  the  same  tip  temperature  as  a  comparable  sized  high- 

TP>c-ic;1"irif~p»     r\  ~rr\V\  e±      it      c  +-  p*  ci  rl  v      ct  qto  ( c  p  p      Pi  nnrp      1/1^  TVi  -i  c 

-^  „  ^^.^v.—.-^w.      j^^^^w      ~~w      o  ~ u  v~  j       .v  . — .. —    ,       ^„w~      -   _  0  i~  _  „      _.y.  i.  ~~^.j 

tip  temperature  can  be  attained  provided  a  proper  value  is 
chosen  for  the  radio- frequency  term  G.   Recall  that  for  the 
rf  probe  the  problem  is  as  follows : 

Differential  equation 


x  '(*'§8 


R2     8R       P°    R"    8t 


Boundary  conditions 
At  R  =  1,  ||  =  0. 


As  R+~,  e+o. 
Initial  condition 
At  t  =  0,  0  =  0. 
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Neglect  blood  flow  and  assume  steady  state 

•Mi)  . 


G_ 

9R       "  R2 


R2  ae  _  g 

R      3R  R    4 

99.  =  G_  +  Cl 

8R  R3    R2 


But  at  R  =  1,  H  =  0  z£>  Cj 


G 


«    • 


89    G     G 
9R   R3    R2 


fl  -    G   +  1    +  r 
u  •    2R2    R   L2 


But  as  R+°°,  6->0  ^=t>C9  =  0 


B  -     G   +  G 

u      2R2    R 


•    • 


T   -  Tn 
At  R  =  1,  8  =  60=  j 


B0  -  2  • 


•    • 


In  Figure  14  0  is  plotted  versus  R.   Once  90  is  chosen, 
one  can  find  a  value  for  G  that  will  result  in  the  steady 
state  temperature  of  the  two  probes  being  the  same  at  the 
probe  tip. 
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c 
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COMPUTER  PROGRAM  1 

PROGRAM  FOR  THE  INSULATED  SPHERE 


100 
150 

2000 


10 
50 


200 


70 

90 

ec 


220 
60 


D I  MENS 
IT  I  ME- 
READ  ( 
FORMAT 
WRITE 
FORMAT 
f 5X,'L 
WRITE 
FORMAT 
TIME  A 
CATION 
EM=1.0 
DELTAU 
DO  10 
TET(K, 
DO  20  C 
TET(  IH 
TET(LI 
X=I 

TET( I , 
)+G*EM 
I  T  I M  E  = 
TAU-IT 
IF( ITI 
IF  (  IT 
ITEST- 
GO  TO 
ITEST= 
GO  TO 
ITEST- 
J  =  0 
DO  6  0 
R=J*DR 

T  F  I   <   T 


ION  TET(500,2) 

0 

5,100) 

(3F 
(6 


1 


(  ' 

1  =  ' 
(6, 

<« 
ND 
(R) 
/(2. 
-EM* 
K  =  l, 
1)=G 

I  =  L 
+  ltl 
-ltl 


10 

50 
i  . 


TA,GtDR,LI 
,2110) 


IH 


13 
00 
', 

AD 
,T 
C  + 
OR 
IH 
.0 
I, 
)  = 
)  = 


BE 
.0 
)B 
T4 

t5 
0) 
T4 

IAL  LOCATIONS' //'  ' ,T3 
^2 

BE 


cTA,DR,G,LI , IH 
,  '6ETA='  ,F4.2,5X 
X,  'IH=« ,14) 


»DR=' ,F6.4,5X, 'G=* , F6.1 


'TEMPERATURE  DISTRIBUTION  AT  SPECIFIED 

,  'TIME(TAU) ' ,T22, ' LO 
'  ) 


'TEMPERATURE(THETA) 
TA*DR**2) 
2 


IH 

TET( IH-1,1) 

TET(LI+1,1) 


2)=EM*((1.-1./X)*TET(I-1,1)+<1.+1./X)*TET(I+1,1) 

*CR**2/<  (X*DR)y<*4) 

ITIME+1 

IME*DELTAU 

ME.GT. 1000)GO  TO  90 

IME.GT. ICO) GO  TO  70 

ITIME/2C*2C 

80 

ITIME/50*5C 

8C 

ITIME/1G0*100 

J=LI ,IH 


T  IMF-  T  TP  <;T  1 


c 


l. in  t 


FORMAT 
TET( J, 
IF  (IT 
GO  TO 
END 


1  1.".  I  T  M  l 

I 


NF.mr-n  Tn 

TCTj  I  2  * 

('  ',F15.4,7XiF6.3',5X,E15.7J 

1  )---TET(  J,  2) 

IME.GT. 4000) STOP 
50 
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c 
c 


COMPUTER  PROGRAM  2 

PROGRAM  FOR  THE  CONSTANT  TEMPERATURE  SPHERE 


100 
150 

2000 


10 

50 

200 


70 

90 

80 


220 
60 


6 
( 

ND 
(R 


D I  MEMS 

ITIME= 

READ  ( 

FORMAT 

WRITE 

FORMAT 

t5X,»L 

WRITE 

FORMAT 

TIME  A 

CATION 

EM=1.C< 

OELTAU 

DO  10 

TET(K, 

TET(  10 

TET( 1C 

DO  20  0 

TET( IH 

X=I 

TET(  I  , 

)+G*EM 

I T  I M  E  = 

TAU=IT 

if ( i  n 

IF  (IT 
ITEST= 

GO  TO 
ITEST= 
GO  TO 
ITEST= 
J  =  0 
LM=LI-1 


ICN  TET(500t2) 

0 

5,10 

(3F 

6,1 


( 


A,G,DR,LI ,IH 

21  10) 

TA,DR,G,LI , IH 

'BETA=«  ,F4.2,5X,  'DR=«  ,  F6.4,5X,  •  G= '  ,F6.1 


,  »IH=»  ,14) 


'TEMPERATURE  DISTRIBUTION  AT  SPECIFIED 


/(2. 
=  EM* 
K  =  l, 
1)  =  0 
fl)  = 
r2)  = 
I  =  L 
+  1,1 


0)BET 

10.0, 

50)  BE 

*,T4, 

13, 5X 

000) 

ST4, 

ADIAL 

,T42, 

C  +  3ET 

DR**2 

IH 

.0 

0.3 

0.3 

I  ,  IH 

)  =  TET( IH-1  ,1) 


LOCATIONS'// '  '  ,T8 
•TEMPERATURE(THETA) 
A*DR**2) 


, 'TIME(TAU) • , 
') 


7  2  2, 'LO 


TO 
TO 


90 
70 


k  =  j  *  D  R 
IF  (  (  I 
WRITE( 
FORMAT 
TET( J, 
IF  (  IT 
GO  TO 
END 


2)=EM*( (  l.-l./X)*TET(I-l,l )  +  ( l.  +  l./X  )*TET  (1  +  1, 

*DR**2/( (X*DR)**4) 

ITIME+1 

IME*OELTAU  • 

ME.GT.100C )G0 

I  ME. GT. 100)60 

ITIME/20*2C 

P r) 

ITIME/50*50 

ITIME/100*100 


T  U 


1  ) 


TIME-ITEST). NE.OJGO  TO  60 
6,220)TAU,R,TET( J, 2) 

('  «,F15.4,7X,F6.3,5XvE15.7) 
1)=TET( J, 2) 
IME.GT.400C) STOP 
50 
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c 
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COMPUTER  PROGRAM  3 

PROGRAM  FOR  THE  INSULATED  CYLINDER 


100 
150 

2000 


10 
50 


200 


70 

90 
80 


2  20 

60 


(3F 
(6,  1 
(  ' 


01  MENS  I CN 
IT  I  ME =0 
READ  (5,10 
FORMAT 
WRITE 
FORMAT 
,5X,'L 
WRITE  (6, 
FORMAT  ( ■ 
TIME  AND 
CATION{R) 
EM=1.0/(2. 
DELTAU=EM* 
DO  10  K=l, 
TET(K,  1)=C 
DO  20  0  I=L 
TET( IH+1,  1 
TET(LI-1,  1 


TET(500,2) 


0  )  3  E 

10.0 

50  )B 

•  ,T4 

13,5. 

000) 

'  »T4 

ADIA 

,T42 

O+BE 

DR** 

IH 

.0 

I,IH 

)=TE 

)=TE 


«DR=» ,F6.4, 5X, «G=« ,F6.1 


TA,G,DR,LI,IH 

,21  10) 

ETA,DR,G,LI ,IH 
,  'BETA='  ,F^fo  2,5X 
IH=« ,14) 

,  'TEMPERATURE 

L  LOCATIONS'//'  '  ,  T3  ,  « T I  ME ( T AU )  •  , T22  , 

, «TEWPERATURE(THETA) ' ) 


DISTRIBUTION 
i 


AT  SPECIFI 


ED 
•LO 


TA: 
2 


DR**2) 


X=I 

TET( 
TET( 


T( IH-1,  1) 
T(LI+1,1 J 


I,2)=E 
1+1,1) 

ITIME=ITIM 
TAU=ITIME* 
IF( ITIME.G 
IF  (ITIME- 
1TEST=ITIM 
GO  TO  80 
ITEST=ITIM 
GO  TO  80 
ITEST=ITIM 
J  =  0 

DO  60  J=LI 
R=J#DR 
I 


M*( ( l.-l./(2. 
)+G*EM/(X**2) 
E+l 

DELTAU 

T. 1000) GO  TO 

GT.100)GO  TO 

E/20*20 

E/50*50 

E/100*100 

,IH 


*X) )*TET( 1-1, 1 )+(!.+!./ (2. *X))* 


90 
70 


l  r  i  V 
WK1 I b( 
FORMAT 
TET( J, 
IF  (IT 
GO  TO 
END 


t  » •  r~ 
1  II  1_' 


t  Tr c  t  \ 

X      I  l.  J  I   I 


MC 


-  0  )  G  Q  T  0  6  0 


(  « 

1  )=T 

I  ME. 
50 


T  I 
I  V 


Ul  I  M  U  t  r\  )  i  i_  i  v  Kt  t  i^ 

',F15.4,7X,F6.3,5X,E15.7) 
E  T  (  J  ,  2  ) 
GTo4000) STOP 
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100 
150 

2000 


10 

50 

200 


70 

90 

80 


2?0 
60 


COMPUTER  PROGRAM  4 
PROGRAM  FOR  THE  CONSTANT 


TEMPERATURE  CYLINDER 


D I M  E  N 
I  TIME 
READ 
FORMA 
WRITE 
FORMA 
,5X,' 
WRITE 
FORMA 
TIME 
CATIO 
EM=1. 
DELTA 
DO  10 
TET(K 
TET(  1 
TET(  1 
DO  20 
TET(  I 
X=I 
TET(  I 
TET(  I 
I T  I M  E 
TAU=I 
IF(  IT 
IF  (  I 
I TE  SI- 
GO  TO 
I  TEST 
GO  TO 
ITEST 
J  =  0 
LM=LI 
DG  60 
R  =  J*U 
IF  (  ( 
WRITE 
FORMA 
TET(  J 
IF  (  I 
GO  TO 
END 


SION 
=  0 

(  5,10 
T  (3F 

(6,1 
T  (  ' 
LI=«f 

(6,2 
TCI 
ANC  R 
N  (  R  )  • 
0/(2. 
U=EM* 

K  =  l, 
il)=0 
0,1)  = 
0,2)  = 
C  I  =  L 
H+l,  1 

,2)=E 
+  1,1) 
=  I T I M 

TIME  * 
I  M  E  .  G 
TIME. 
=  ITIM 

80 
=  IT  I  M 

SO 
=  ITIM 

-1 

O-Li'i 

R 

I  T  I  M  E 
(6,22 
T  (  ' 
,  1  )=T 

T  I  ME. 
50 


TET(50C,2) 


C)BETA,G,DP,LI ,IH 
21  10) 
TA,DR,G 
■BETA=« 
,  •  IH=' 


1  u  »  u  f 
50  )  B  E 
•  ,T4, 
13, 5X 
000) 
'  ,T4, 
ADIAL 
,T42 

0  +  BET 
DR**2 
IH 

.0 

1.0 

1.0 

1  ,IH 

)  =  TET 


,  L I  ,  I  H 
,F4.2,5X, 


14) 


DR=',F6.4,5X,  'G=«  ,F6.1 


•TEMPERATURE  DISTRIBUTION  AT  SPECIFIED 

LCC AT  IONS' //'  »  ,T8,  « TIME(TAU) « ,T22,  » LO 
•  TEMPERATURE (THETA)  ' ) 
A#DR**2) 


( IH-1,1) 


M*((l-.-l./(2.*X) 

)+G*EM/( X**2) 

E+l 

DELTAU 

T.  1000)  GO  TO  90 

GT.IC'OIGO    TO    70 

E/20*20 

E/50*50 

E/100*1C0 


)*TET( I-1,1)+(1.+1./(2.*X)>* 


T    I    I 

f  A  (1 


60 


-ITEST) c NE.0JGC1    TO 

0)TAU,P,TET( J, 2) 

'  ,F) 5.4, 7X,F6.3,5X,E15 

E  T  (  J  ,  2  ) 

GT.4C0U  STOP 


7) 


90 


c 

c 


COMPUTER  PROGRAM  5 

PROGRAM  FOR  THE  HIGH-RESISTANCE  SPHERE 


100 

150 

2000 

50 

12 

11 


220 
200 


13 
14 


DIMENS 

READ  ( 

FORMAT 

WRITE 

FORMAT 

tI5) 

FORMAT 

TIME  A 

CATION 

SQBET- 

TAU=0„ 

SQTAU== 

DR=0.2 

BTAU=B 

SQBTAU 

00  200 

R=I*DR 

EX=SQ3 

EXNEG^= 

A=(R-1 

B=(R-1 

TET(I) 

RFC( B) 

WRITE 

FORMAT 

CO NT  IN 

IF{TAU 

GO  TO 

TAU=TA 

TAU=TA 

ITAU=T 

IF(  ITA 

GO  TO 


ION 


p,100) 


TET( 500) 

BETA,THETO,  IH 
(2F10. 0,110) 

(6,150)  BETA,THETO,IH 
{  ■  «  ,T4,  'BETA='  ,  F8.4,5X,  •THETO=«  ,  F  10.  1 ,  5X,  •  I  H=» 


(  •  1 ' ,  T4 t  ' 

ND    RAOIAL 
(R>  »  ,T42,  ' 
SQRT(BETA) 
2 
SQRT(TAU) 


TEMPERATURE  DISTRIBUTION 
LOCATIONS1//*  »»T8 

) 


, .ot'TIME(TAU)« ,T22, 
TEMPERATURE(THETA) • 


AT    SPECIFIED 
•LO 


cmp\ 


ETA*TAU 
=SQRT(BTAU) 
1  =  1,  IH 

ET*(R-l.O) 

-EX 

.0 )/( 2. 0*SOTAU)- SQBTAU 

•0)/(2.0*SQTAU)+SQBTAU 

=THETO*( 1,/(2.J*R)*<EXP< EXN EG ) *ERFC ( A )  +  EXP (  EX ) *E 

)  ) 

(6,220)    TAUr.R,TET(  I  ) 

{'     ' ,F15.4, 7X,F6.3,5X,E15.7) 
UE 

cGT.2.0)     GO    TO    13 
14 

U+0.5 
U  +  0*2 
AU 

U.GT.30)     STOP 
50 
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TABLE  I . 

Values  for  6/G  for  various  values  of  non-dimensional 
time  (t)  and  non-dimensional  location  (R) :   insulated 
sphere . 

3  =  0.0 
t  0.1     0.5     1.0     2.0     5.0     10.0    20.0    120.0  * 


R 


1.0  0.0501  0.1320  0.1825  0.2375  0.3085  0.3554  0.3940  0.4943 

1.2  0.0408  0.1190  0.1685  0.2228  0.2933  0.3401  0.3786  0.4748 

1.4  0.0280  0.0974  0.1441  0.1967  0.2661  0.3125  0.3508  0.4620 

1.6  0.0180  0.0761  0.1192  0.1693  0.2369  0.2827  0.3208  0.4340 

1.8  0.0114  0.0580  0.0968  0.1438  0.2093  0.2543  0.2921  0.4051 

2.0  0.0074  0.0434  0.0776  0.1213  0.1842  0.2283  0.2651  0.3650 

2.2  0.0049  0.0322  0.0617  0.1018  0.1618  0.2050  0.2420  0.3500 

2.4  0.0034  0.0237  0.0487  0.0851  0.1421  0.1842  0.2207  0.3252 

2.6  0.0024  0.0174  0.0385  0.0701  0.1248  0.1657  0.2016  0.3050 

2.8  0.0018  0.0128  0.0300  0.0589  0.1096  0.1492  0.1846  0.2852 

•7      r\  r\      rvn-iT  r>      nr>or  r>      n  i  v  *  n      n/ion  o      r\  r\  c  ->  n      i   7  a  c  r\      -\  n  c\  ■?  r>      icon 

O.U  U.UUX.J  O.UUC.J  \J   •  <J  t*  *J -t  u.Utu^  O.Uw'O^  u.J-^/-rvJ  O.jlvj^o  u  .   i^Oi-u 

3.4  0.0008  0.0054  0.0143  0.0334  0.0742  0.1097  0.1431  0.2450 

3.8  0.0005  0.0032  0.0087  0.0227  0.0571  0.0898  0.1216  0.2225 

4.2  0.0003  0.0020  0.0054  0.0154  0.0439  0.0736  0.1039  0.2050 

4.6  0.0002  0.0013  0.0034  0.0104  0.0336  0.0604  0.0891  0.1900 

5.0  0.0002  0.0009  0.0023  0.0071  0.0258  0.0496  0.0766  0.1720 

7.0  0.0000  0.0002  0.0005  0.0012  0.0065  0.0184  0.0369  0.1254 

10.0  0.0000  0.0001  0.0001  0.0002  0.0009  0.0039  0.0126  0.0980 

15.0  0.0000  0.0000  0.0000  0.0001  0.0002  0.0006  0.0039  0.0540 


*   Steady  State  Value 
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R 


TABLE  II. 

Values  for  6/G  for  various  values  of  non-dimensional 
time  (t)  and  non-dimensional  location  (R) :   insulated 
sphere . 

3  =  0.1 
t  0.1     0.5     1.0     2.0     5.0     10.0    15.0    20.0 


1.0  0.0499  0.1296  0.1765  0.2241  0.2756  0.2988  0.3058  0.3085 

1.2  0.0406  0.1168  0.1627  0.2097  0.2608  0.2840  0.2909  0.2936 

1.4  0.0278  0.0953  0.1388  0.1843  0.2346  0.2575  0.2645  0.2671 

1.6  0.0179  0.0744  0.1145  0.1578  0.2068  0.2294  0.2363  0.2390 

1.8  0.0113  0.0566  0.0926  0.1333  0.1807  0.2030  0.2098  0.2124 

2.0  0.0073  0.0423  0.0741  0.1118  0.1573  0.1791  0.1859  0.1885 

2.2  0.0049  0.0313  0.0587  0.0934  0.1368  0.1581  0.1648  0.1673 

2.4  0.0034  0.0231  0.0463  0.0777  0.1189  0.1397  0.1462  0.1488 

2.6  0.0024  0.0170  0.0363  0.0645  0.1054  0.1235  0.1300  0.1325 

2.8  0.0018  0.0125  0.0284  0.0534  0.0899  0.1094  0.1157  0.1182 

3.0  0.0013  0.0092  0.0222  0.0441  0.0782  0.0970  0.1032  0.1057 

3.4  0.0008  0.0052  0.0135  0.0300  0.0592  0.0767  0.0826  0.0850 

3.8  0.0005  0.0031  0.0082  0.0201  0.0449  0.0609  0.0665  0.0688 

4.2  0.0003  0.0020  0.0137  0.0137  0.0340  0.0485  0.0538  0.0560 

4.6  0.0002  0.0013  0.0092  0.0092  0.0258  0.05S8  0.0438  0.0459 

5.0  0.0002  0.0009  0.0063  0.0063  0.0195  0.0311  0.0357  0.0377 

7.0  0.0000  0.0002  0.0011  0.0011  0.0048  0.0105  0.0134  0.0149 

10.0  0.0000  0.0001  0.0002  0.0002  0.0007  0.0021  0.0033  0.0039 

15.0  0.0000  0.0000  0.0000  0.0000  0.0001  0.0003  0.0007  0.0010 
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TABLE  III. 


Values  for  9/G  for  various  values  of  non-dimensional 
time  (t)  and  non-dimensional  location  (R) :   insulated 
sphere . 


3  =  1.0 


t  0.1  0.5  1.0  2.0  5.0  10.0  * 
R 

1.0  0.0479  0.1108  0.1354  0.1439  0.1527  0.1528 

1.2  0.0389  0.0989  0.1230  0.1363  0.1400  0.1401 

1.4  0.0266  0.0795  0.1023  0.1152  0.1189  0.1190 

1.6  0.0171  0.0611  0.0821  0.0943  0.0978  0.0980 

1.8  0.0108  0.0458  0.0647  0.0762  0.0796  0.0797 

2.0  0.0070  0.0339  0.0504  0.0611  0.0644  0.0644 

2.2  0.0046  0.0249  0.0391  0.0488  0.0519  0.0520 

2.4  0.0032  0.0182  0.0302  0.0393  0.0419  0.0419 

2.6  0.0023  0.0133  0.0233  0.0311  0.0339  0.0339 

2.8  0.0017  0.0098  0.0180  0.0249  0.0274  0.0275 

3.0  0.0013  0.0072  0.0138  0.0199  0.0222  0.0223 

3.4  0.0007  0.0041  0.0085  0.0128  0.0148  0.0148 
3.8  0.0005  0.0024  0.0050  0.0083  0.0099  0.0100 

4.0  0.0004  0.0019  0.0040  0.0067  0.0082  0.0082 

4.5  0.0002  0.0011  0.0023  0.0040  0.0051  0.0051 
5.0  0.0002  0.0007  0.0014  0.0024  0.0032  0.0033 

7.0  0.0000  0.0001  0.0003  0.0005  0.0007  0.0007 

10.0  0.0000  0.0000  0.0001  0.0001  0.0001  0.0001 

15.0  0.0000  0.0000  0.0000  0.0000  0.0000  0.0000 


Steady  state  value 
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TABLE  IV. 


Values  for  9/G  for  various  values  of  non-dimensional 
time  (t)  and  non-dimensional  location  (R) :   insulated 


sphere . 


3  =  10.0 

t  0.1     0.5     1.0  * 

R 

1.0  0.0337  0.0433  0.0434 
1.2  0.0265  0.0356  0.0357 
1.4  0.0175  0.0253  0.0254 
1.6  0.0110  0.0172  0.0173 
1.8  0.0069  0.0116  0.0116 

•2.0  0.0045  0.0078  0.0079 
2.2  0.0030  0.0054  0.0054 
2.4  0.0021  0.0038  0.0038 
2.6  0.0015  0.0027  0.0027 
2.8  0.0011  0.0019  0.0019 


.5  .  U   U.UUUO   U.UU-LH   U.Uult 

4.0  0.0003  0.0004  0.0004 
7.0  0.0000  0.0000  0.0000 


Steady  state  value 
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TABLE  V 


Values  for  6/90  for  various  values  of  non-dimensional 
time  (x)  and  non-dimensional  location  (R) ;   constant 
temperature  sphere. 


3  =  0.0 
t0.2     0.6     1.0     2.0     5.0     10.6    30.9 


R 


1.0  1.0000  1.0000  1.0000  1.0000  1.0000  1.0000  1.0000 

1.2  0.6265  0.7126  0.7396  0.7670  0.7913  0.8045  0.8164 

1.4  0.3765  0.5107  0.5552  0.6011  0.6424  0.6648  0.6853 

1.6  0.2142  0.3649  0.4196  0.4776  0.5309  0.5602  0.5870 

1.8  0.1144  0.2584  0.3176  0.3829  0.4446  0.4789  0.5105 

2.0  0.0569  0.1806  0.2397  0.30S5  0.3759  0.4140  0.4494 

2.2  0.0263  0.1242  0.1801  0.2493  0.3202  0.3612  0.3994 

2.4  0.0119  0.0839  0.1342  0.2016  0.2741  0.3171  0.3578 

2.6  0.0044  0.0554  0.0992  0.1629  0.2357  0.2800  0.3226 

2.8  0.0016  0.0358  0.0725  0.1314  0.2033  0.2485  0.2925 

3.0  0.0005  0.0226  0.0524  0.1057  0.1757  0.2213  0.2664 

3.4  0.0000  0.0084  0.0170  0.0677  0.1215  0.1771  0.2236 

4.0  0.0000  0.0015  0.0044  0.0334  0.0766  0.1287  0.1757 

5.0  0.0000  0.0001  0.0003  0.0091  0.0345  0.0770  0.1222 

8.0  0.0000  0.0000  0.0000  0.0001  0.0021  0.0161  0.0467 

20.0  0.0000  0.0000  0.0000  0.0000  0.0000  0.0000  0.0008 
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TABLE  VI 


Values  for  6/8Q  for  various  values  of  non-dimensional 
time  (x)  and  non-dimensional  location  (R) :   constant 
temperature  sphere. 


3  =  0.1 
t0.2      0.6      1.0      2.0      5.0      10.6     30.9 


1, 

,0 

1, 

,0000 

1, 

,0000 

1, 

,0000 

1, 

.0000 

1, 

,0000 

1, 

,0000 

1, 

,0000 

1. 

,2 

0. 

,6238 

0, 

,7069 

0, 

,7319 

0 

.7557 

0, 

,7735 

0, 

,7799 

0, 

,7822 

1, 

,4 

0, 

,3736 

0, 

.5032 

0, 

,5444 

0 

.5842 

0, 

,6144 

0, 

,6254 

0, 

,6293 

1. 

,6 

0, 

.2121 

0, 

.3575 

0, 

,4081 

0 

.4585 

0, 

.4974 

0, 

,5117 

0, 

,5168 

1, 

,8 

0, 

,1130 

0, 

,2520 

0, 

,3066 

0 

.  3633 

0, 

,4084 

0, 

,4251 

0, 

,4331 

2, 

,0 

0, 

.0562 

0, 

,1754 

0, 

,2300 

0 

.2898 

0, 

,3389 

0, 

,  3575 

0, 

.3642 

2, 

.2 

0, 

,0259 

0, 

.1202 

0, 

,1718 

0 

.2319 

0, 

.2835 

0, 

,3035 

0. 

,3107 

2, 

,4 

0, 

,0110 

0 

.0809 

0, 

,1275 

0 

.1859 

0 

.2387 

0, 

,2597 

0, 

.2673 

.2, 

.6 

0, 

,004  3 

0, 

.0533 

0, 

,0938 

0 

.1490 

0, 

.2019 

0, 

,2236 

0, 

.2316 

2, 

,8 

0, 

,0015 

0, 

.0344 

0, 

,0683 

0 

.1193 

0 

.1715 

0, 

.1936 

0, 

,2018 

3, 

.0 

0, 

,0005 

0, 

.0217 

0, 

.1)492 

0 

.0953 

0  , 

,  1461 

0  , 

,16  85 

,  JL  /  U  / 

3, 

.4 

0, 

,0000 

0 

.0080 

0 

.0246 

0 

.0601 

0 

.1066 

0, 

,1287 

0 

,1373 

4, 

.0 

0. 

,0000 

0 

.0015 

0, 

.0079 

0 

.0293 

0 

.0670 

0 

.0877 

0 

.0964 

5, 

.0 

0, 

,0000 

0 

.0000 

0 

.0009 

0 

.0078 

0 

.0307 

0 

.0478 

0 

.0561 

8, 

.0 

0 

.0000 

0 

.0000 

0 

.0000 

0 

.0000 

0 

.0023 

0 

.0081 

0 

.0153 

20, 

,0 

0, 

,0000 

0 

.0000 

0 

.0000 

0 

.0000 

0 

.0000 

0 

.0000 

0, 

.0001 
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TABLE  VII 


Values  for  6/60  for  various  values  of  non-dimensional 
time  (t)  and  non-dimensional  location  (R) :   constant 
temperature  sphere. 


3  =  1.0 
t  0.2      0.6      1.0      2.0      5.0      10.6  * 


R 


1 . 

,0 

1 

.0000 

1, 

.0000 

1 

.0000 

1, 

,0000 

1 

.0000 

1, 

.0000 

1, 

,2 

0 

.6002 

0, 

.6614 

0 

.6740 

0, 

.6809 

0 

.6822 

0, 

,6823 

1. 

,4 

0 

.3491 

0 

.4442 

0 

.4648 

0, 

.4764 

0 

.4788 

0, 

,4788 

1. 

,6 

0, 

.1938 

0, 

.3000 

0 

.3253 

0, 

.3399 

0, 

.3429 

0, 

,3431 

1. 

,8 

0, 

.1015 

0, 

.2024 

0 

.2297 

0, 

.2461 

0 

.2496 

0, 

.2496 

2. 

,0 

0 

.0498 

0 

.1356 

0 

.1629 

0, 

.1801 

0 

.1839 

0 

.1839 

2. 

,2 

0 

.0227  • 

0 

.0899 

0 

.1156 

0, 

.1329 

0 

.1368 

0 

.1369 

2, 

,4 

0 

.0096 

0 

.0588 

0 

.0820 

0 

.0987 

0 

.1027 

0 

.1027 

2, 

,6 

0 

.0037 

0, 

.0378 

0 

.0579 

0, 

,0  736 

0 

.0776 

0 

.0777 

2, 

,8 

0, 

.0013 

0 

.0239 

0 

.0406 

0 

.0551 

0 

.0590 

0, 

.0590 

3 , 

,  n 

0 

,  n  0  0  4 

n 

.0148 

0 

.0283 

0 

.0413 

0 

.0451 

0 

.0451 

4 . 

,  0 

0 

*~*  *"*  r\  ^ 

.  uuuu 

V 

r\  o  r\  o 

,UUU3 

A 

on  .1  n 

n 

p  no  q 

.  U  V  *J   w 

n 

ni  o a 

>  \^   -«-  >~-    i 

q 

.  0124 

8, 

.0 

0 

.0000 

0 

.0000 

0 

.0000 

0 

.0000 

0 

.0001 

0 

.0001 

Steady  state  value 
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TABLE  VIII. 


Values  for  0/Go  for  various  values  of  non-dimensional 
time  (t)  and  non-dimensional  location  (R) :   constant 
temperature  sphere. 


3   =  10.0 

t  0.2      0.6      1.0 
R 


1 , 

,0 

1.0000 

1.0000 

1. 

.0000 

1. 

,2 

0.4384 

0.4427 

0, 

.4427 

1. 

,4 

0.1951 

0.2016 

0 

.2016 

1, 

,6 

0.0866 

0.0937 

0 

.0937 

1, 

,8 

0.0378 

0.0442 

0, 

.0443 

2, 

,0 

0.0160 

0.0211 

0 

.0212 

2, 

,4 

0.0025 

0.0049 

0 

.0050 

3. 

,0 

0.0000 

0.0006 

0 

.0006 

4, 

,0 

0.0000 

0.0000 

0 

.0000 

!•    r>  j .1,. 

~  *-    n   4-   ~    Vn1 
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TABLE  IX. 

Values  for  0/G  for  various  values  of  non-dimensional 
time  (t)  and  non-dimensional  location  (R) :   insulated 
cylinder . 

3  =  0.0 

t  0.1     0.5     1.0     2.0  5.0  10.0  15.0  20.0 
R 

1.0  0.0694  0.2449  0.3985  0.6241  1.0637  1.5190  1.8401  2.0936 

1.2  0.0628  0.2342  0.3862  0.6104  1.0475  1.5029  1.8257  2.0770 

1.4  0.0517  0.2123  0.3597  0.5780  1.0130  1.4662  1.7861  2.0388 

1.6  0.0411  0.1867  0.3274  0.5415  0.9682  1.4181  1.7365  1.9883 

1.8  0.0327  0.1664  0.2935  0.4998  0.9182  1.3656  1.6800  1.9307 

2.0  0.0263  0.1381  0.2606  0.4576  0.8660  1.3058  1.6198  1.8692 

2.2  0.0216  0.1176  0.2298  0.4165  0.8.134  1.2468  1.5580  1.8057 

2.4  0.0180  0.1000  0.2017  0.3774  0.7615  1.1877  1.4957  1.7416 

2.6  0.0152  0.0854  0.1767  0.3409  0.7112  1.1295  1.4338  1.6776 

2.8  0.0131  0.0731  0.1543  0.3072  0.6629  1.0724  1.3729  1.6144 


3.0  0.0113  0.0631  0.1356  0.2764  0.6168  1.0171  1.3133  1.5524 
3.4  0.0C33  0.0480  C. 10-17  0.2232  0.5320  0.9122  1.1990  1.4326 
3.8  0.0070  0.0376  0.0821  0.18U4  U.4S/1  0.8157  1.0921  1.3196 


4.2  0.0037  0.0301  0.0655  0.1465  0.3919  0.7276  0.9930  1.2137 

4.6  0.004S  0.0249  0.0532  0.1198  0.3356  0.6480  0.9015  1.1151 

5.0  0.0040  0.0209  0.0442  0.0990  0.2874  0.5762  0.8175  1.0255 

7.0  0.0020  0.0104  0.0213  0.0453  0.1372  0.3186  0  .  4964" 0  .  6616 

10.0  0.0010  0.0051  0.0102  0.0209  0.0576  0.1598  0.2404  0.3496 

15.0  0.0005  0.0024  0.0049  0.0104  0.0290  0.0688  0.1233  0.1945 
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TABLE  X. 

Values  for  0/G  for  various  values  of  non-dimensional 

time  (i)  and  non-dimensional  location  (R) :   insulated 
cylinder . 

3  =  0.1 

t  0.1     0.5  1.0  2.0  5.0  10.0  15.0  20.0 

R 

1.0  0.0690  0.2396  0.3823  0.5771  0.8913  1.1131  1.2072  1.2521 

1.2  0.0625  0.2290  0.3702  0.5637  0.8769  1.0983  1.1922  1.2371 

1.4  0.0514  0.2074  0.3444  0.5345  0.8446  1.0650  1.1586  1.2034 

1.6  0.0409  0.1823  0.3129  0.4977  0.8033  1.0219  1.1152  1.1598 

1.8  0.0325  0.1575  0.2801  0.4582  0.7576  0.9741  1.0667  1.1111 

2.0  0.0262  0.1347  0.2483  0.4184  0.7105  0.9242  1.0162  1.0603 

2.2  0.0214  0.1146  0.2187  0.3798  0.6636  0.8741  0.9652  1.0091 

2.4  0.0179  0.0975  0.1919  0.3434  0.6179  0.8248  0.9150  0.9585 

2.6  0.0152  0.0832  0.1679  0.3095  0.5739  0.7769  0.8660  0.9092 

2.8  0.0130  0.0713  0.1469  0.2784  0.5322  0.7308  0.8188  0.8616 

3.0  0.0113  0.0615  0.1287  0.2500  0.4928  0.6868  0.7735  0.8159 

3.4  0.0088  0.0-167  0.0993  0.2014  0.4212  0.6053  0.6892  0.7306 

3.8  0.0070  0.0366  0.0778  0.1625  0.3591  0.5325  0.6133  0.6536 


4.2  0.0057  0.0295  0.0621  0.1318  0.3059  0.4679  0.5455  0.5846 

4.6  0.0047  0.0243  0.0506  0.1078  0.2605  0.4111  0.4851  0.5229 

5.0  0.0040  0.0204  0.0420  0.0891  0.2222  0.3611  0.4316  0.4680 

7.0  0.0020  0.0102  0.0203  0.0409  0.1054  0.1918  0.2435  0.2726 

10.0  0.0010  0.0049  0.0097  0.0189  0.0447  0.0835  0.1124  0.1316 

15.0  0.0004  0.0023  0.0047  0.0094  0.0225  0.0413  0.0568  0.0690 
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TABLE  XI . 

Values  for  6/G  for  various  values  of  non-dimensional 
time  (t)  and  non-dimensional  location  (R) :   insulated 
cylinder . 

6  =  1.0 
t   0.1      0.5      1.0      2.0      5.0      10.0  * 


R 


1.0  0.0661  0.1988  0.2729  0.3269  0.3480  0.3487 

1.2  0.0598  0.1892  0.2626  0.5162  0.3372  0.3379 

1.4  0.0491  0.1701  0.2412  0.2938  0.3147  0.3154 

1.6  0.0390  0.1484  0.2162  0.2673  0.2877  0.2884 

1.8  0.0330  0.1273  0.1909  0.2400  0.2601  0.2607 

2.0  0.0249  0.1083  0.1672  0.2140  0.2334  0.2341 

2.2  0.0204  0.0919  0.1456  0.1899  0.2087  0.2094 

2.4  0.0171  0.0780  0.1266  0.1682  0.1863  0.1869 

2.6  0.0145  0.0664  0.1100  0.1488  0.1661  0.1668 

2.8  0.0124  0.0569  0.0958  0.1316  0.1482  0.1488 

3.0  0.0108  0.0490  0.0836  0.1166  0.1323  0.1329 

3.4  0.0084  0.0373  0.0643  0.0919  0.1060  0.1065 

5.8  0.0067  0.C295  0.0504  0.0732  0.0856  0.0861 

4.2  0.0054  0.0236  0.0405  0.0590  0.0698  0.0705 

4.6  0.0045  0.0195  0.0529  0.0482  0.0575  0.0580 

5.0  0.005S  0.0165  0.0274  0.0400  0.0480  0.0484 

7.0  0.0019  0.0082  0.0154  0.0819  0.0226  0.0228 

10.0  0.0010  0.0040  0.0064  0.0089  0.0104  0.0105 

15.0  0.0004  0.0019  0.0051  0.0045  0.0051  0.0052 


Steady  state  value 
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TABLE  XII. 

Values  for  6/G  for  various  values  of  non-dimensional 
time  (t)  and  non-dimensional  location  (R) :   insulated 
cylinder . 

3  =  10.0 

x   0.1      0.5  1.0  * 

R 

1.0     0.0449    0.0633  0.0656 

1.2     0.0400    0.0578  0.0581 

1.4     0.0323    0.0487  0.0490 

1.6     0.0255    0.0399  0.0402 

1.8     0.0208    0.0326  0.0328 

2.0     0.0163   '0.0267  0.0269 

2.2     0.0133    0.0221  0.0222 

2.4     0.0111    0.0184  0.0186 

2.6     0.0094    0.0156  0.0158 

2.8     0.0081    0.0133  0.0135 

3.0     0.0071    0.0116  0.0117 

3.4     0.0055    0.0089  0.0C90 

3.8     0.0044    0.0071  0.0071 

4.2     0.0036    0.0058  0.0058 

4.6     0.0030    0.0048  0.0048 

5.0     0.0025    0.0040  0.0041 

7.0     0.0013    0.0020  0.0021 

10.0     0.0006    0.0010  0.0010 

15.0     0.0005    0.0005  0.0005 

*   Steady  state  value 
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TABLE  XIII. 

Values  for  8/60  for  various  values  of  non-dimensional 
time  (t)  and  non-dimensional  location  (R) :   constant 
temperature  cylinder. 

$  =  0.0 

t  0.1     0.5     1.0     2.0     5.0     10.0    15.0  20.0 

R 

1.0  1.0000  1.0000  1.0000  1.0000  1.0000  1.0000  1.0000  1.0000 

1.2  0.6091  0.7771  0.8214  0.8543  0.8856  0.9027  0.9108  0.9159 

1.4  0.3266  0.5939  0.6724  0.7319  0.7890  0.8205  0.8355  0.8449 

1.6  0.1510  0.4448  0.5469  0.6273  0.7058  0.7494  0.7703  0.7834 

1.8  0.0590  0.3254  0.4412  0.5368  0.6329  0.6870  0.7129  0.7292 

2.0  0.0190  0.2319  0.3524  0.4583  0.5683  0.6313  0.6617  0.6809 

2.2  0.0050  0.1608  0.2783  0.3898  0.5107  0.5813  0.6156  0.6372 

2.4  0.0010  0.1083  0.2171  0.3302  0.4590  0.5360  0.5737  0.5976 

2.6  0.0001  0.0707  0.1672  0.2784  0.4123  0.4947  0.5354  0.5612 

2.8  0.0000  0.0447  0.1270  0.2534  0.3702  0.4569  0.5001  0.5277 

3.0  0.0000  0.0273  0.0951  0.1946  0.3320  0.4221  0.4676  0.4967 

3.4  0.0000  0.0092  0.0509  0.1327  0.2659  0.3604  C.4093  0.4412 

7      Q     n     nn  r\n     r\     on  o  n     n     o  ?  c  ^     n     nooi      o     911  c     n     7iT7^     n     7.  r  o  o  n     -7  n  n  a 

•J    •     VJ  \J     .     \J   \J    \J   \J  W     i    VU    <-    /  \J     •     O    tt   kJ  \J  »k/«Ww*->-i-  v      .      (_     J .     j      v '  \J     •     xj  KJ     f     *J  \J     •     %^    ^j    W   \J  \J     •     ^/  ^/    *_,   ~\ 

4.2  0.0000  0.0007  0.0120  0.0569  0.1669  0.2618  0.3145  0.3495 

4.6  0.0000  0.0001  0.0053  0.0356  0.1306  0.2224  0.2755  0.3114 

5.0  0.0000  0.0000  0.0022  0.0216  0.1011  0.1883  0.2411  0.2774 

7.0  0.0000  0.0000  0.0000  0.0011  0.0239  0.0770  0.1200  0.1534 

10.0  0.0000  0.0000  0.0000  0.0000  0.0015  0.0157  0.0368  0.0586 

15.0  0.0000  0.0000  0.0000  0.0000  0.0000  0.0010  0.0066  0.0175 
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TABLE  XIV. 

Values  for  8/60  for  various  values  of  non-dimensional 
time  (x)  and  non-dimensional  location  (R) :   constant 
temperature  cylinder. 

3  =  0.1 

t  0.1     0.5     1.0  2.0  5.0  10.0  15.0  20.0 

R 

1.0  1.0000  1.0000  1.0000  1.0000  1.0000  1.0000  1.0000  1.0000 

1.2  0.6072  0.7713  0.8126  0.8412  0.8669  0.8725  0.8749  0.8758 

1.4  0.3249  0.5859  0.6590  0.7106  0.7576  0.7679  0.7724  0.7740 

1.6  0.1500  0.4365  0.5316  0.6013  0.6659  0.6803  0.6864  0.6888 

1.8  0.0585  0.3180  0.4258  0.5087  0.5877  0.6056  0.6132  0.6162 

2.0  0.0188  0.2259  0.3379  0.4297  0.5203  0.5411  0.5501  0.5536 

2.2  0.0049  0.1561  0.2654  0.3620  0.4617  0.4850  0.4952  0.4990 

2.4  0.0010  0.1048  0.2060  0.3040  0.4103  0.4350  0.4470  0.4512 

2.6  0.0001  0.0683  0.1580  0.2542  0.3650  0.3923  0.4044  0.4090 

2.8  0.0000  0.0431  0.1195  0.2115  0.3249  0.3537  0.3665  0.3714 

3.0  0.0000  0.0262  0.0892  0.1752  0.2893  0.3193  0.3328  0.3380 

5.4  0.0000  0.0088  0.0475  0.1180  0.2293  C.2G09  0.2754  0.2811 

3.8  0.0000  0.0026  0.0238  0.0776  0.1814  0.2137  0.2289  0.2349 

5.0  0.0000  0.0000  0.0020  0.0186  0.0876  0.1176  0.1332  0.1397 

7.0  0.0000  0.0000  0.0000  0.0009  0.0228  0.0420  0.0546  0.0606 

15.0  0.0000  0.0000  0.0000  0.0000  0.0000  0.0004  0.0020  0.0038 
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TABLE  XV. 


Values  for  6/90  f°r  various  values  of  non-dimensional 
time  (t)  and  non-dimensional  location  (R) :   constant 
temperature  cylinder. 

3   =  1.0 


R 


0.1 


0.5 


1.0 


2.0 


Steady  state  value 


5.0 


1.0 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.2 

0.5908 

0.7246 

0.7466 

0.7549 

0.7566 

1.4 

0.3102 

0.5220 

0.5607 

0.5757 

0.5789 

1.6 

0.1414 

0.3716 

0.4221 

0.4422 

0.4466 

1.8 

0.0547 

0.2604 

0.3175 

0.3414 

0.3467 

2.0 

0.0175 

0.1790 

0.2381 

0.2645 

0.2706 

2.2 

0.0045 

0.1203 

0.1778 

0.2054 

0.2121 

2.4 

0.0009 

0.0789 

0.1319 

0.1598 

0.1668 

2.6 

0.0001 

0.0504 

0.0970 

0  .1243 

0.1316 

2.8 

0.0000 

0.0  312 

0.0708 

0.0967 

0.1040 

3.0 

0.0000 

0.0188 

0  .0511 

0.0752 

0.0825 

7  A 

•J     a  '  I 

0.0000 

0  .  0  0  6 1 

0.02  57 

U  .  0  4  5  2 

0.0521 

3.8 

0.0000 

0.0018 

0.0123 

0.0269 

0.0331 

5.0 

0.0000 

0.0000 

0.0009 

0.0051 

0.0086 

7.0 

0.0000 

0.0000 

0.0000 

0.0002 

0.0009 

15.0 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 
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TABLE  XVI. 

Values  for  6/60  for  various  values  of  non-dimensional 
time  (t)  and  non-dimensional  location  (R) :   constant 
temperature  cylinder. 

3  =  10.0 

t   0.1       0.5       1.0  * 
R 


1.0 

1.0000 

1.0000 

1.0000 

1.2 

0.4631 

0.4889 

0.4890 

1.4 

0.2028 

0.2421 

0.2421 

1.6 

0.0809 

0.1208 

0.1210 

1.8 

0.0284 

0.0607 

0.0609 

2.0 

0.0085 

0 .0306 

0.0308 

2.2 

0.0021 

0.0155 

0.0157 

2.4 

0.0004 

0.0078 

0.0080 

2.6 

0.0000 

0.0039 

0.0041 

2.8 

0.0000 

0.0020 

0.0021 

3.0 

0.0000 

0.0010 

0  ,  f)f)  l  0 

3.4 

0.0000 

0.0002 

0.0003 

3.8 

0.0000 

0.0001 

0.0001 

Steady  state  value 
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